Critical Casimir effect in classical binary liquid mixtures 
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If a fluctuating medium is confined, the ensuing perturbation of its fluctuation spectrum gener- 
ates Casimir-like effective forces acting on its confining surfaces. Near a continuous phase transition 
of such a medium the corresponding order parameter fluctuations occur on all length scales and 
therefore close to the critical point this effect acquires a universal character, i.e., to a large extent 
it is independent of the microscopic details of the actual system. Accordingly it can be calculated 
theoretically by studying suitable representative model systems. We report on the direct measure- 
ment of critical Casimir forces by total internal reflection microscopy (TIRM), with femto-Newton 
resolution. The corresponding potentials are determined for individual colloidal particles fioating 
above a substrate under the action of the critical thermal noise in the solvent medium, constituted 
by a binary liquid mixture of water and 2,6-lutidine near its lower consolute point. Depending 
on the relative adsorption preferences of the colloid and substrate surfaces with respect to the two 
components of the binary liquid mixture, we observe that, upon approaching the critical point of the 
solvent, attractive or repulsive forces emerge and supersede those prevailing away from it. Based on 
the knowledge of the critical Casimir forces acting in film geometries within the Ising universality 
class and with equal or opposing boundary conditions, we provide the corresponding theoretical 
predictions for the sphere — planar wall geometry of the experiment. The experimental data for 
the effective potential can be interpreted consistently in terms of these predictions and a remarkable 
quantitative agreement is observed. 

PACS numbers: 05.70. Jk,82.70.Dd,68.35.Rh 



I. INTRODUCTION 



Fluctuation-induced forces 



At macroscopic scales thermal or quantum fluctuations 
of a physical property of a system are typically negligible 
because fluctuations average out to zero upon increasing 
the length and time scales at which the system is stud- 
ied. At the micro- and nano-meter scale, instead, fluctu- 
ations become generally relevant and, if externally con- 
trolled and spatially confined, they give rise to novel phe- 
nomena. An example thereof is provided by the Casimir 
force acting on conducting bodies [l| , which is due to the 
confinement of quantum fluctuations of the electromag- 
netic field in vacuum and which influences the behavior of 
micrometer-sized systems ranging from colloids to micro- 
and nano-electromechanical systems (MEMS, NEMS). 

Thermal fluctuations in condensed matter typically oc- 
cur on a molecular scale. However, upon approaching the 
critical point of a second-order phase transition the fluc- 
tuations of the order parameter cf) of the phase transition 
become relevant and detectable at a much larger length 
scale ^ and their confinement results in a fluctuation- 
induced Casimir force /c acting on the confining sur- 
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faces This so-called critical Casimir force /c has a 
range which is set by the correlation length ^ of the fluc- 
tuations of the order parameter. Since near the critical 
point ^ can reach up to macroscopic values, the range 
of /c can be controlled and varied to a large extent by 
minute temperature changes close to the critical temper- 
ature Tc. We shall show that this control of the ther- 
modynamic state of the system is a manageable task. 
This implies that the critical Casimir force can be easily 
switched on and off, which allows one to identify it rela- 
tive to the omnipresent background forces. In addition, 
by proper surface treatments of the confining surfaces, 
the force can be relatively easily turned from attractive 
to repulsive [1, Q in contrast to the Casimir force stem- 
ming from electromagnetic fluctuations for which such a 
change requires carefully chosen hulk materials provid- 
ing the solid walls and the fluid in between 0. Such 
a repulsive force might be exploited to prevent stiction 
in MEMS and NEMS, which would open signiflcant per- 
spectives for applications. Finally, at Tc the strength 
of the critical Casimir force can easily compete with or 
even dominate dispersion forces, with which it shows the 
same algebraic decay, however without suffering from the 
weakening due to retardation effects. The universality 
of /c means that the same force is generated near the 
critical point of liquid-vapor coexistence of any fluid or 
near the consolute point of phase segregation of any bi- 
nary or multicomponent liquid mixture. This allows one 
to pick and use those representatives of the universality 
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class which in addition optimize desired performances of 
MEMS and NEMS. This provides a highly welcome flex- 
ibility. 

The fluctuation-induced forces generated by confining 
the fluctuations of electromagnetic fields in the quantum 
vacuum (Casimir effect) or of the order parameter in a 
critical medium (critical Casimir effect) have a common 
description within the field-theoretical approach. Ac- 
cordingly, the connection between these two effects goes 
well beyond the mere analogy and it indeed becomes an 
exact mapping in some specific cases of spatial dimen- 
sion d, geometries, and boundary conditions. This deep 
connection actually justifies the use of the term "criti- 
cal Casimir force" when referring to the effective force 
due to the confinement of critical fluctuations. On the 
other hand, from a theoretical point of view the quan- 
tum and the critical Casimir effect are also distinct in 
that the quantum one in vacuum corresponds to a free 
field theory whereas the critical one is described by a 
more challenging non-Gaussian field theory. 



B. Finite-size scaling 

The theory of finite-size scaling (see, e.g., Refs. [1, 01) 
predicts that in the vicinity of the critical Casimir 
force /c and its dependence on temperature are described 
by a universal scaling function which depends only on 
the gross features of the system and of the confining 
surfaces, i.e., on the so-called universality class of the 
phase transition occurring in the bulk and on the ge- 
ometry and surface universality classes of the confining 
surfaces d, 0, Hi- The latter characterize the bound- 
ary conditions (BC) [1, H, 0, the surfaces impose on 
the fluctuations of the order parameter of the underlying 
second-order phase transition. The actual physical na- 
ture of the order parameter (j) depends on which kind of 
continuous phase transition is approached: in the case we 
shall be mainly concerned with in the following, i.e., the 
consolute point of phase segregation in binary liquid mix- 
tures, (j) is given by the difference between the local and 
the mean concentration of one of the two components of 
the mixture (see, c.f.. Sec. If 11 B] for further details). For 
binary liquid mixtures the confining surfaces generically 
exhibit preferential adsorption of one of the two compo- 
nents of the mixture, resulting in an enhancement of the 
order parameter (j) close to the surface. (This amounts 
to the presence of symmetry-breaking surface fields, see, 
e.g., Refs. d, 0, 111-) One usually refers to the corre- 
sponding boundary conditions as (-I-) or (— ) depending 
on whether the surface favors > or < 0, respec- 
tively. Due to its universal nature, the critical Casimir 
force can be studied via representative models which are 
amenable to theoretical investigations. Since due to uni- 
versality microscopic details can only in a rather limited 
way be blamed for potential discrepancies, the resulting 
predictions face very stringent experimental tests. 

Most of the available theoretical and experimental 



studies focus on the film geometry in which the system 
undergoing the second-order phase transition is confined 
between two parallel surfaces of large transverse area S 
at a distance L. For this geometry and assuming that 
the only relevant thermodynamic variable is the temper- 
ature T (possible additional variables such as the concen- 
tration are set to their critical values), renormalization- 
group theory shows 9, f 0] that the critical Casimir force 
/c scales as 

= ^.mL/iifh (1) 

in three spatial dimensions (d = 3), where dix) is a uni- 
versal scaling function, and r is the reduced deviation 
from the critical temperature Tc such that t > corre- 
sponds to the disordered (homogeneous) phase. If, as it 
is usually the case, the homogeneous phase is located at 
high temperatures in the phase diagram of the system, 
one defines r = (T — Tc)/Tc. However, there are also 
cases - such as the one we shall be interested in (see, 
c.f.. Fig. [8]) - in which this phase is located at low tem- 
peratures so that there one defines r = — (T — Tc)/Tc. 
The system-specific (i.e., non- universal) amplitudes 
in Eq. ([T]) enter into the algebraic behavior of the bulk 
correlation length ^ of the order parameter (j) upon ap- 
proaching the critical point: 

ar-^0^)=^^\rr . (2) 

In what follows we shall mainly consider = Co"; 
which forms with S,q a universal amplitude ratio U^^^^ = 
/£,Q ~ f.9 [ni, Ell in those cases in which ^(t < 0) 
is finite. (Renormalization-group theory tells that in the 
bulk there are only two independent non-universal am- 
plitudes, say ^J" and = (0)/(— t)^ of the order pa- 
rameter below Tc, all other non- universal amplitudes can 
be expressed in terms of them and universal amplitude 
ratios [ll]. Here (3 is the critical exponent which charac- 
terizes the singular behavior (0) ~ (—t)^ of the average 
order parameter {(j)) for r — > 0^, with /? = 0.3265(3) 
for the three-dimensional Ising universality class ^11|.) 
The bulk correlation length ^ can be inferred from, e.g., 
the exponential decay of the two-point correlation func- 
tion of the order parameter. The algebraic increase of 
^ [Eq. ^] is characterized by the universal exponent v 
which equals 0.630f(4) for the three-dimensional Ising 
universality class [lli] which captures, among others, the 
critical behavior of binary liquid mixtures close to the 
demixing point as studied experimentally here. 

C. Theoretical predictions and previous 
experiments 

For the Ising universality class with symmetry- 
breaking boundary conditions theoretical predictions for 
the universal scaling function d are available from field- 
theoretical [lE[i3 and Monte Carlo studies [isl. [isl. . 



3 



The critical Casimir force turns out to be attractive for 
equal boundary conditions (BC) on the two surfaces, i.e., 
(+,+) or (—,—), whereas it is repulsive and generically 
stronger for opposing boundary conditions, i.e., (+, — ) 
or (— , +). In the presence of such boundary conditions, 
for topographically or chemically 18 1 patterned con- 
fining surfaces or for curved surfaces 1^ [20| theoretical 



results are available primarily within mean-field theory. 

Following theoretical predictions and suggestions [2ll |. 
previous indirect evidences for both attractive and repul- 
sive critical Casimir force were based on studying fluids 
close to critical endpoints (see Ref. [22| for a more de- 
tailed summary). Under such circumstances, the film 
geometry with parallel planar walls can be indeed ex- 
perimentally realized by forming complete wetting fluid 
films [2^ in which a liquid phase is confined between 
a solid substrate (or another spectator phase) and the 
interface with the vapor phase and its thickness L can 
be tuned by undersaturation, in particular off critical- 
ity. Upon changing pressure and temperature one can 
drive the liquid film towards a second-order phase tran- 
sition which nonetheless keeps the confining liquid-vapor 
interface sharp. The fluctuations of the associated order 
parameter, confined within the film of thickness L, give 
rise to a critical Casimir pressure (related to i? [2l|) which 
acts on the liquid-vapor interface, displacing it from the 
equilibrium position it would have under the effect of 
dispersion forces alone, i.e., in the absence of critical 
fluctuations. This results in a temperature-dependent 
change of L. Based on the knowledge of the relation- 
ship between L and pressure, by monitoring this vari- 
ation it is possible to infer indirectly the magnitude 
of the Casimir force which drives this change of thick- 
ness. This approach has been used for the study of 
wetting films of *He at the normal-superfluid transi- 
tion 241, for '^He-'*He mixtures close to the tricritical 
point [25|], and for classical binary liquid mixtures close 
to demixing transitions [1^ [13]. The film thickness L 
has been determined by using capacitance 0, HI] or X- 
ray reflectivity measurements l2al . or ellipsometry [27l |. 
For the results of Refs. [25||. and [26] the quan- 

titative agreement with the theoretical predictions for 
the correspon ding bulk and surface universality classes 
(see Refs. H Efl p, p, [M H, S [M!, [H, HI], 
|33 . l34ll . and [IJ, [la, [la, [33 . r espectively) are excel- 
lent |24j or remarkably good [25l. [2a ] . For ''He [13] one 
has Dirichlct-Dirichlet boundary conditions, for '^He-'*He 
mixtures ||25j Dirichlct-(+) boundary conditions, and in 
Ref. [2^ (-I-, — ) boundary conditions hold. 



D. Direct determination of critical Casimir forces 

The aim of the experimental investigation discussed 
here is to provide a direct determination of the Casimir 
force, by measuring the associated potential ^c- On di- 
mensional grounds and on the basis of Eq. ([1]), the scale 
of this potential is set by k^Tc and therefore, as real- 



ized in Ref. [26[, in order to enhance the strength of 
the critical Casimir force it is desirable to engage crit- 
ical points with higher Tc compared to those of the A- 
transition investigated in Refs. [13, [l^- This consider- 
ation suggests classical fluids as natural candidates for 
the critical medium. The experimentally driven prefer- 
ence for having Tc and the critical pressure to be close to 
ambient conditions can be satisfied by numerous binary 
liquid mixtures which exhibit consolute points for phase 
segregation. From Eq. ([1]) one can infer a rough estimate 
of the critical Casimir force fc- For an object which ex- 
poses an effective area S ^ 1 /im^ to a wall at a distance 
L = 100 nm, and for = 300 K one finds fc < 4pN. 
Since the scaling function t?(x) vanishes upon moving 
away from criticality, i.e., ^{\x\ — > 00) — > 0, and because 
one is interested in probing also larger distances L, one 
needs force measurements with a force resolution which 
is significantly better than pN. Atomic force microscopy 
at room temperature cannot deliver fN accuracy. This 
required sensitivity can, however, be achieved by using 
total internal reflection microscopy (TIRM), which en- 
ables one to determine the potential of the effective forces 
acting on a colloidal particle near a wall, by monitoring 
its Brownian motion in a solvent. Choosing as the sol- 
vent a suitable binary liquid mixture allows one to in- 
vestigate the critical Casimir force on the particle which 
arises upon approaching the demixing transition of the 
mixture. Such a second-order phase transition falls into 
the bulk universality class of the Ising model. In this 
geometrical setting the fluctuation spectrum of the criti- 
cal medium (i.e., the binary liquid mixture) is perturbed 
by the confinement due to a flat wall and by the pres- 
ence of the spherical cavity. The curvature of one of the 
two confining surfaces introduces an additional length 
scale and thus leads to an extension of the scaling form 
in Eq. ([T]) such that the scaling function t9 additionally 
depends on the ratio between the radius R of the colloid 
and the minimal distance z between the surface of the col- 
loid and the flat surface of the substrate Jc.f., Sec. Ill Al 
here 2 plays the role of L in Eq. ([T])]. At present, for 
arbitrary values of z and radii of curvature, theoretical 
predictions for the critical Casimir force in a geometri- 
cal setting involving one non-planar surface are available 
only within mean-field theory, both for spherical 19, 2^] 
and ellipsoidal [33] particles, which demonstrate that the 
results of the so-called Derjaguin approximation are valid 
for z/R < 1 [H, [10] (see, c.f.. Sec. [III. Beyond mean- 
field theory and for various universality classes, theoret- 
ical results have been obtained in the so-called protein 
limit corresponding to z/i?, ^/z ;2> 1 [1, [3^] where R in- 
dicates the typical size of the, in general nonspherical, 
particle. However, at present this protein limit is not 
accessible by TIRM because for small particles far away 
from the substrate (through which the evanescent optical 
field enters into the sample) the signal of the scattered 
light from the particle is too weak. The experimentally 
relevant case is the opposite one of a large colloidal par- 
ticle close to the wall. Although in d = 3 theoretical 
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results for the full scaling function of the sphere-plate 
geometry are not available, in this latter case one can 
take advantage of the Derjaguin approximation in order 
to express the critical Casimir force Fq acting on the col- 
loid in terms of the force acting within a film geometry, 
which was investigated successfully via Monte Carlo sim- 
ulations in Refs. [15|, [l^. This is explained in detail in 
Sec. Ill Al in which we present the theoretical predictions 
for the scaling function of the critical Casimir force (and 
of the associated potential) for the case of a sphere near a 
wall immersed into a binary liquid mixture at its critical 
composition. On the other hand, in Sec. IIIBl we discuss 
the expected behavior of the effective potential of the 
colloid if the binary liquid mixture is not at its critical 
concentration so that, upon changing the temperature, 
it undergoes a first-order phase transition. The discus- 
sions in Sec. [iTl form the basis for the interpretation of 
the experimental results. The experimental setting is de- 
scribed in Sec. IIIII In Sec. IIII Al we recall the principles 
of TIRM and of the data analysis, whereas in Sec. IIII Bl 
we discuss the specific choice of the binary mixture used 
here and how one can experimentally realize the various 
boundary conditions. In Sec. II VI we present in detail the 
experimental results, comparing them with the theoret- 
ical predictions, for mixtures both at critical and non- 
critical compositions. A summary and a discussion of 
perspectives and of possible applications of our findings 
are provided in Sec. |Vl Part of the analysis presented 
here has been reported briefly in Ref. [s^]. (Fo r a peda- 
gogical introduction to the subject see Ref. |38j.) 

II. THEORETICAL PREDICTIONS 
A. Critical composition 

1. General properties 

The critical Casimir force Fq acting on a spherical par- 
ticle of radius R, at a distance z of closest approach 
from the flat surface of a substrate and immersed in 
a near-critical medium at temperature T ~ Tc takes, 
for stron g p referential adsorption, the universal scaling 
form [ll,l0,[3i 

The scaling function K^^'^\x, A) depends, in addition, 
on the combination of (sphere, plate) [(s,p)] boundary 
conditions imposed by the surfaces of the sphere and of 
the plate and on the phase from which the critical point is 
approached (i.e., on the sign of r, with K^'^^ correspond- 
ing to r ^ 0). (In line with Eq. ^ and with the standard 
notation in the literature, the one-phase region is denoted 
by -|- and the two-phase region by — . These signs should 
not be confused with the signs (-I-, +) etc. indicating, also 
in line with the literature, the character of the boundary 



conditions of the two confining surfaces (s,p). In order 
to avoid a clumsy notation we suppress or use these two 
notations in a selfevident way.) The scaling form of the 
associated potential ^c{z) = /^""ds Fc{s) follows by in- 
tegration of Eq. ([3]). In the two limiting cases A 1 
and A -C 1 it is possible to calculate K{x, A) on the ba- 
sis of the so-called small-sphere expansion and Derjaguin 
approximation, respectively [19L |20, 39']. In the former 
case one finds in three space dimensions, for r > and 
symmetry-breaking boundary conditions (s,p) — (±,-l-) 
(see Eq. (7) in Ref. which also includes higher-order 
terms) 

+ ^ ' ^ c+ 2/3/'^ +^ ' (4) 

+ 0(A-2/3/^-i), 

where P/v ~ 0.518. In this limit, the force acting 
on the "small" particle is determined, to leading order, 
by the interaction between the particle and the aver- 
age order parameter profile {4>{z))oo/2 induced by the 
planar wall in the absence of the particle, i.e., in a 
semi- infinite system (00/2). This profile is character- 
ized for T > by the universal scaling function P+ en- 
tering Eq. (g]): {(j){z))oa/2 = {'t>)oo-r<QP+{z/C)i whcrc 
{4')oo,-T<o — C^T^ is the value of the order parameter in 
the bulk (00) corresponding to the reduced temperature 
—T — > 0~. The universal constant c+ in Eq. ^ char- 
acterizes the critical adsorption profile P+{x ^ 0) — > 
c+x~^/^ , whereas a = A^^/Bcj, is the universal ratio [3§| 
between the non-universal amplitudes A^f, and of the 
critical order parameter profile in the semi-infinite system 
{(t>{z))oc/2.T=Q = A^{2z)~^^^ and of the two-point corre- 
lation function in the bulk ((^(r)</)(0))oo,r=o = -B^r^^'^/", 
respectively. In turn, Arf, (and therefore B^) can be ex- 
pressed in terms of the two independent non-universal 
amphtudes ^(J" and via A^ = c+C^/{2^'^)-P/'' . (For 
a detailed discussion of the values of these universal am- 
plitude ratios we refer the reader to Refs. [191. i20.. .39.].) 

2. Derjaguin approximation 

Equation @ is useful to discuss the behavior of col- 
loids which arc small compared to their distance from the 
plate. However, in the experiment discussed in Sec. IIIII 
the distance z is typically much smaller than the radius 
R of the particle. This case can be conveniently dis- 
cussed within the Derjaguin approximation, which yields 
in three dimensions 19] 

K{x,A-.Q) ^ A-'^d{x) , (5) 

where the expression for 'd{x) is determined further below 
in terms of the scaling function d of the critical Casimir 
force /c acting within a film [see Eq. ([1])]. 

The scaling functions _)(a::) and i?(+__|_)(a;) for the 
boundary conditions (-1-,—) and (-t-,-l-) relevant to the 
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FIG. 1: (a) Geometry of the Derjaguin approximation for the 
plate-sphere geometry, (b) cross section through the center 
of the sphere and normal to the plate. In (a), the base of 
a cap of the sphere with radius r{9) is shown as a thin line. 
The inner circle of the grey ring dS{9) is this thin line shifted 
by RsinOdO towards the center of the sphere. In (b), 5 = 
r{9 + dO) — r{8) is the cross section of the grey ring dS{6) 
shown in (a). 



study of the critical properties of binary liquid mixtures 
at their critical compositions have been determined by 
Monte Carlo simulations [l^[l^. Within the Derjaguin 
approximation, valid for A <C 1, i.e., if the radius R of 
the colloid is much larger than the minimal separation z 
between the surface of the colloid and the flat substrate, 
the curved surface of the colloid is considered to be made 
up of successive circular rings of inflnitesimal area dS{9) 
and radius r{9) which are parallel to the substrate and 
are at a normal distance L{9) = z + R{1 — cos 6') from an 
opposing identical circular ring on the surface of the sub- 
strate (see Fig. [1} . Assuming additivity, the contribution 
dFc of each single pair of rings to the total Casimir force 
Fc is given by 



dFc _ dS{0) 



(6) 



where ■(? is the scaling function of the critical Casimir 
force acting within the film geometry [see Eq. ([ij]. Here 
it is convenient to express d not as a function of u = 
t(L/^o)^^'^ as in Eq. ([T|) but as a function of a; = L/^ 
where x — u'^ for u > and x = Ue (—uY for u < 
with J7{g„p ~ 1.9 [111 for the three-dimensional Ising uni- 
versality class we are interested in. The radius r(9) of 
the ring is given by r{9) = Rsm9 and therefore its area 
is dS{9) = n[r{9 + d9)]^ - T:[r{9)]^ = 2ttR^ sin9 cos9 d9. 
The total force Fq is obtained by summing all the contri- 
butions dFc{9) of the circular rings, up to the maximal 
angle 9m 



fcRT 



0=0 



dS{9) 
LH9) 



(7) 



= ""72 is a natural choice, neglecting any influences 
from the back side of the sphere. However, we shall see 



below that its specific value does not affect the result 
in the limit R ^ z. For R ^ z the integral (due to 
the denominator) is dominated by the contributions it 
picks up at small angle 9, so that one can approximate 
L{9)/z ~ 1 + {R/z){9^/2) and therefore 



Fciz) 27ri?2 



knT 



d9 



[1 + ^9^2] 



rm + -0V2]z/O. 



[For 6'm — 00 this is identical with Eq. (4) in Ref. [19].] 
Introducing the variable Z = 1 -t- {R/ z){9'^ /2), one can 
write the previous expression as 



Fc{z) 2ttR 



2 , dlj^^lz/0 (9) 



where Im = I + (R/ z){9l^/2). In the limit R/ z oo, 
— > cxD independently of 9m , so that the integral can be 
extended up to oo and 



where 



d{x) = 27r / dl -j^^ilx 



(10) 



(11) 



The potential ^c{z) associated with the Casimir force is 
given by 



= -27: r dv (\ - \\ ^{v z/0 (12) 



kBT 



= - e{z/0 

z 

where we have changed the variable I i-^ v = 
ly, exchanged the order of the remaining integrals 
Jj^ dy Jy dv = /j^ dv dy, and introduced the scaling 
function 

e(x) = 2^^ dv(^^ - ^^^{vx). (13) 

According to Eqs. ^TU\i and p2p . for separations z much 
smaller than the radius of the colloid, the Casimir force 
and the Casimir potential increase linearly upon increas- 
ing the radius R of the colloid. At the bulk critical point, 
,9(0) = 9(0) = 7n9(0) and d'{0) = 27r??'(0), whereas 
0'(O) = oo. If in the film geometry the force is at- 
tractive (repulsive) at all temperatures, within the Der- 
jaguin approximation the same sign holds also in the 
sphere-plate geometry. At the critical concentration the 
Casimir force acting on a (-I-) colloid in front of a (— ) 
substrate is the same as the one acting on a (— ) colloid 
in front of a (-I-) substrate. This is no longer true for 
non-critical concentrations. Although the Derjaguin ap- 
proximation is expected to be valid only for i? 3> z, the 
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comparison between the results of the mean-field calcu- 
lation [l^, for the actual sphere-plate geometry and 
the ones of the corresponding Derjaguin approximation 
based on the mean-field theory (MFT) for the film geom- 
etry show good agreement even for z/ R up to 0.4 0.5. 



3. Theoretical predictions for scaling functions 

For the universality class of the three-dimensional Ising 
model, the scaling functions for the Casimir force in the 
film geometry - which enter into Eq. ~ have been 
determined in Refs. [H, [3 for {+, +) and {+, -) BC 
(or, equivalently, (— , — ) and (— , BC) by Monte Carlo 
simulations. Due to the presence of strong corrections 
to scaling, the amplitudes of the corresponding numer- 
ical estimates for i?(+__|_)(a;) and are affected 
by a systematic uncertainty of about 20% [H, [l^ . The 
numerical data presented in Refs. [IB, [lB| are very well 
fitted by certain analytic ansatze (at least in the range 
of scaling variable which has been investigated numeri- 
cally) which, in turn, can be used in order to calculate 
the corresponding scaling functions and _ i 
for the potential (Fig. [H see also Fig. 2(d) in Ref. [33|) 
as well as and for the force (Fig. [3|). The 
simulation data for the film scaling functions '(?(_|_ +)(a;) 
and _)(a;) can actually be fitted even by functions 
of various shapes (the asymptotic behavior of which for 
large |x| is, however, fixed, see further below). This leads 
to different estimates of the scaling functions outside the 
range of the scaling variable for which the Monte Carlo 
data are currently available. This results also in different 
estimates of and obtained 
from (x) and (x) via Eqs. and PT|) . 
However, the uncertainty of the estimates for the shapes 
is negligible compared to the inherent systematic uncer- 
tainty associated with the amplitudes of i?(+.+)(a;) and 
-){x). For a detailed discussion of these issues we 
refer to Ref. 

The critical Casimir force fc{T,L) between two pla- 
nar walls [see Eq. ([1])] with symmetry-breaking bound- 
ary conditions is expected to vary as exp(— L/$) as a 
function of i ^ for r > (see, e.g., Ref. [401 ^^'^ 
particular the footnote Ref. [23] therein). Accordingly, 
^{+,±){x » 1) = A± cc^e"^ and from Eqs. ^ and ^ 
one finds 

^{+,:lz){x > 1) = 2ttA± e"^ and ^^^^ 
0(+,±)(a; > 1) = 27rA±xe"'', 

for the critical Casimir force and potential, respectively, 
in the sphere-plate geometry. The analysis of the Monte 
Carlo data presented in Figs. 9 and 10 of Ref. [l^ yields 
A^^ = -1.51(2) and A^^ ^ 1.82(2), respectively, for the 
data sets therein indicated as {i) whereas it yields a'"^^ = 

— 1.16(2) and A^^'^ = 1.38(2) for the corresponding data 
set (ii). (We recall here that the data sets (i) and (m) 



turn out to be proportional to each other, see Refs. [15|, 
[3 for details.) 

Figure [3] shows that the critical Casimir force for the 
sphere-plate geometry exhibits the same qualitative fea- 
tures as in the film geometry: For (-I-, -|-) [(-I-, — )] BC the 
force is attractive (repulsive) and attains its maximum 
strength for r > (r < 0), corresponding to the one- 
phase (two-phase) region. For fixed values of the scaling 
variable, the strength of the repulsive force for (+, — ) BC 
is larger than the one of the attractive force in the case 
of (-I-, -|-) BC. The inset of Fig. [3] compares the estimate 
for the scaling function _|_)(a; = L/^) — up to its nor- 
malization '(^(U-4-)(0) — based on the Monte Carlo data 
of Refs. (Til llq (solid line) with the early estimate of 
Ref. [l^, which is based on the pointwise and linear in- 
terpolation between the exactly known film scaling func- 
tions in d — 2 and d = 4, such as to obtain an estimate of 
i?(+.+)(a; — L/S) for d = 3 (dashed line). Although this 
latter estimate captures correctly some qualitative fea- 
tures of the universal scaling function i?(_|__+)(x = i/^), 
it fails to be quantitatively accurate, as the comparison 
with the Monte Carlo estimate reveals. The same consid- 
eration applies to the corresponding estimates for 

Equations and p3|) . together with Fig.[21 form the 
theoretical basis for the interpretation of the experimen- 
tal results for the effective interaction potential between 
a spherical colloidal particle and a planar wall, immersed 
into a binary liquid mixture at its critical composition 
and near its consolute point. 



4. Deviations from strong adsorption 

The theoretical analyses presented above and in 
Refs. [H, [H, [131 assume that the confining surfaces are 
characterized by a sufficiently strong preferential adsorp- 
tion for one of the two components of the mixture, cor- 
responding to (-I-) or (— ) fixed-point boundary condi- 
tions in the sense of renormalization-group theory [^ [3] . 
Within the coarse-grained field-theoretical description of 
the binary mixture close to a boundary B in terms of the 
order parameter [1,[3|, the preferential adsorption is ac- 
counted for by a surface contribution —hg Jg dS 0(x G B) 
to the effective free energy of the system, where the "sur- 
face field" hg summarily quantifies the strength of the 
preferential adsorption. Indeed, /is > [hg < 0] fa- 
vors (/> > [0 < 0] at the boundary B so that, for 
\hs\ large enough, \(t){z)\ oc z~^/'^ at normal distances 
z ^ (but still large on molecular scales) from B 
The (-I-) and (— ) boundary conditions correspond to the 
limits hs — > -l-oo and —00, respectively, of strong pref- 
erential adsorption. Within this coarse-grained descrip- 
tion the gross features of the relation between hg and the 
material properties of the wall and the mixture can be 
inferred from the behavior of experimentally accessible 
quantities such as critical adsorption profiles or excess 
adsorption (see, e.g., Refs. [4l|, |42|). For a weak adsorp- 
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FIG. 2: Scaling functions 0(+__f,) and 0(+ _) of the Casimir 
potential $c [see Eq. dHJ] for (+,+) and (+, -) BC, re- 
spectively, within the Derjaguin approximation and for the 
three-dimensional Ising universality class, as functions oi u — 
t{z/^q)^^'^ with u ~ 0.630. The thick solid and dashed lines 
have been obtained via Eq. (I13|l on the basis of the Monte 
Carlo estimates for (+ _) presented in Refs. p^ . [T^ . 

indicated by (i) and (ii), respectively, in Figs. 9 and 10 of 
Ref. (The thick soUd Unes agree with the estimates re- 

ported in Fig. 2(d) of Ref. 113].) attains its minimum 

value S(™'"'j — —3.6 (solid line) and —2.8 (dashed line) both 
for iimin — 0.54, whereas attains (smoothly) its max- 

imum value — 19 (solid line) and 14 (dashed line) 

both for Umin — —0.03. The first derivatives of 0(+,_) and 
0(+,+) diverge logarithmically for m ^ 0. The thin lines for 
w > 1 indicate the asymptotic behaviors of 0(u ^ 1) given 
in Eq. (|14[) with the numerical values of the coefficients A± 
indicated from top to bottom for the corresponding curves. 
For (-I-, -I-) boundary conditions the asymptotic expressions 
are indistinguishable from 0(+.+)(ti) for it > 4. 



tion preference, the corresponding hs might be so small 
that upon approaching the critical point one effectively 
observes a crossover in the kind of boundary condition 
imposed on the order parameter. The critical Casimir 
force reflects such [43] or related [i^l crossover behav- 
iors; in the film geometry, depending on the film thick- 
ness, the force can even change sign [H, On the 
basis of scaling arguments one expects that for moder- 
ate adsorption preferences the scaling function in Eq. ([T]) 
additionally depends on the dimensionless scaling vari- 
ables Us^i = Uihs^iL^^/'^ , i = 1,2, where /is,i, /is, 2 are 
the effective surface fields at the two confining surfaces, 
Oi > are corresponding non-universal constants, and 
Ai ~ 0.46 is the so-called surface crossover exponent at 
the so-called ordinary surface transition 0, . One can 
associate a length scale ii = {ai\hs^i\)~'^^^^ with each 
surface field, such that the theoretical predictions dis- 
cussed before are valid for L ^ ii, i.e., ys^i ±oo, 
whereas corrections depending on ti/L are expected to 
be relevant for L ~ ii. For ii 3> L, instead, the prefer- 
ential adsorption of the wall i is so weak that a crossover 
occurs towards boundary conditions which preserve the 



-5 t , ^ , , , I 

-10 -5 5 10 15 20 

FIG. 3; Scaling functions '?{+,+) and '?(+,-) of the Casimir 
force Fc [see Eq. (fTOjl ] for (-I-, -I-) and (-I-, — ) BC, respec- 
tively, within the Derjaguin approximation and for the three- 
dimensional Ising universality class, as functions oi u = 
Tiz/^+Y^" with u ~ 0.630. The thick solid and dashed lines 
have been obtained via Eq. ((TTJ on the basis of the Monte 
Carlo estimates for presented in Refs. (T^, [l^ . 

indicated by (i) and {ii), respectively, in Figs. 9 and 10 of 
Ref. [l6l |. '?(-!-,+) attains its minimum value i^'™'"''^ ~ —4.9 
(solid line) and —3.8 (dashed line) both for Wmin — 2.6, 
whereas attains its maximum value ~ 21 (solid 

line) and 16 (dashed line) both for iimin — —1.5. The second 
derivatives of i9(+_+)(m) and _)(ii) diverge logarithmically 
for u — > 0. The thin lines for u > 1 indicate the asymptotic 
behaviors "diu ^ 1) given in Eq. p4[) with the numerical 
values of the coefficients A± indicated from top to bottom 
for the corresponding curves. For (+, +) boundary condi- 
tions the asymptotic expressions are indistinguishable from 
i9(+,+)(m) for u > 5. In the inset we compare the estimate 
for (a;)/i?(_|. .(.)(0) as a function of a; = z/^ based on the 

Monte Carlo data of Refs. [H, (solid fine, MC) with the 
one presented in Ref. and obtained by interpolating lin- 
early and pointwise the exactly known film scaling functions 
in d = 2 and d = 4 in order to obtain an estimate for d = 3 
(dotted line, Interpol.). The Monte Carlo estimate for this 
ratio is the same for both data sets (i) and {ii) in Ref. [l^. 

I— > — symmetry and there appears to be no effective 
enhancement of the order parameter upon approaching 
the wall. Heuristically, the length scales ii can be in- 
terpreted as extrapolation lengths Zcx,i oc ii in the sense 
that for small enough li ^ Q the order parameter profile 
behaves as 0)| ~ (z z^^^i)-!^/" [1, IH, S^j upon 

approaching the wall i. Within the concept of an extrap- 
olation length the effects of a physical wall with a moder- 
ate preferential adsorption (which implies ii ^ 0) on the 
order parameter are equivalent to those of a fictitious wall 
with strong preferential adsorption (which means ii = 0) 
displaced by a distance — Zox,i from the physical wall. Al- 
though t his p icture is consistent only within mean-field 
theory [1, it turns out to be useful for the interpreta- 
tion of experimental results [4l| and simulation data [i^j 
as an effective means to take into account corrections to 
the leading critical behavior. Assuming that this carries 
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over to the critical Casimir forces, a film of thickness L 
and moderate adsorption at the confining surfaces is ex- 
pected to be equivalent to a film with strong adsorption 
and thickness Zox.i + i + 2^ox,2 > L. On the same footing, 
a sphere of radius R and a plate at a surface-to-surface 
distance z, both with moderate preferential adsorption, 
should behave as a sphere of smaller radius R — Zox.sph 
and a plate at a distance Zex,pi + ^ + ^ox,sph > z, both with 
strong preferential adsorption. We anticipate here that 
the interpretation of the experimental data presented in 
Sec. IIVBI does not require to account for the effect de- 
scribed above, even though we cannot exclude the pos- 
sibility that such corrections might become detectable 
upon comparison with theoretical data with a smaller 
systematic uncertainty than the ones considered here. 

B. Noncritical composition 

1. General properties 

In this section we consider thermodynamic paths ap- 
proaching the critical point from the one-phase region 
by varying the temperature at fixed off- critical compo- 
sitions, e.g., ca ^ c^, where ca is the concentration 
of the A component of a mixture. For systems with 
a lower consolute point these paths lie below the up- 
wards bent phase boundary of first-order phase transi- 
tions in the temperature-composition (T, c) parameter 
space (see, c.f., the vertical paths in Fig. [SJb)). Per- 
forming experiments along such paths is another use- 
ful and interesting probe of the critical Casimir force, 
because the corresponding Casimir scaling function ac- 
quires an additional scaling variable S — sgn{h)L/lh, 
where Ih = Io\{ca — c^a) I c-'a\~^ ^ ^ ^^^id Iq is a nonuniversal 
amplitude. The bulk field h is proportional to the differ- 
ence (^[lA — Ms) ~ (ma — fJ'B)c of the chemical potentials 
of the two components of a binary liquid mixture. If this 
difference is nonzero one has ca 7^ c\ for species A in 
the bulk. The nonuniversal amplitude can be deter- 
mined from the corresponding correlation length Ih which 
is experimentally accessible by measuring the scattering 
structure factor for various concentrations ca > c\ at 
T = Tc- This nonuniversal amplitude is actually related 
to the two independent nonuniversal amplitudes and 
(7^ (see the discussion below Eq. (HI) by the expression 

where S and 7 are standard bulk critical exponents and 
Q2 and R-^ are universal amplitude ratios [ll|, [l^, l48j 
leading to [Q2/{6R^)Y/^ « 0.38 in d = 3. 

So far, for the sphere-plane geometry of the present ex- 
periment there are no theoretical results available for the 
critical Casimir force for thermodynamic states which lie 
off the bulk critical composition. However, based on the 
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FIG. 4: Schematic phase diagram of a binary liquid mixture 
with a lower demixing transition point in terms of tempera- 
ture T and concentration ca of the A species. The solid curve 
encloses the two-phase region separating via first-order phase 
transitions the a and /? phases rich in A and B species, re- 
spectively, terminating at the critical point CP. The dashed 
line indicates the first-order bridging phase transition which 
occurs if a fluid mixture is conflned between a planar wall and 
a sphere of radius R and at distance L possessing the same 
adsorption preference, here for the a phase. The bridging 
transition ends at the critical point Q and separates a region 
in the bulk phase diagram in which a phase preferable by 
walls condenses and forms a bridge connecting the wall and 
the sphere, from the region in the bulk phase diagram where 
such a bridge is absent. Although the bridging transition is 
a (quasi-) first-order phase transition, in the bulk phase dia- 
gram it is described by a line instead of a coexistence region, 
because it is an interfacial phase transition. 

theoretical analysis of the critical Casimir force for films 
[13, mi and sphere-sphere geometries 20], we expect that 
along suitably chosen paths of fixed off-critical composi- 
tions the critical Casimir force is strongly influenced by 
capillary bridging transitions. Moreover, if the bulk field 
h is nonzero (-I-, -t-) and (— , — ) BC are no longer equiva- 
lent. 



2. Bridging transition 

A bridging transition is the analogue of capillary con- 
densation [50j for geometries in which one or both sur- 
faces are non-planar. (However, there is a conceptual dif- 
ference. Whereas capillary condensation corresponds to 
an actual shift of the bulk phase diagram, bridging tran- 
sitions are interfacial phase transitions which leave the 
bulk phase diagram unchanged but can be described as 
if effectively the bulk phase boundary of first-order phase 
transitions is shifted [51.,^].) It occurs at temperatures 
for which two phases may exist, i.e., for T above Tc in 
the case of a binary liquid mixture with a lower conso- 
lute point, and it depends on the adsorption properties 
of the surfaces. If, say, both surfaces favor the a phase 
rich in species A over the (3 phase rich in species B, one 
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expects the a phase to form a bridge between the sur- 
faces for some chemical potential \xa of species A such 
that [I A < IJ-'a^ where /x™ the value corresponding to 
bulk coexistence. Alternatively, this occurs at a concen- 
tration (mole fraction) ca < Ca slightly smaller than its 
value Ca at bulk coexistence. If the surfaces favor the 
(3 phase, the (3 phase fills the gap between the surfaces 
forming a bridge for fj.A > y^X^ phase separation 

line for this morphological transition occurs on the other 
side of the bulk phase diagram, i.e., for ca > c™ (Fig. [5]). 

Bridging may occur in the presence of thin wetting lay- 
ers on both surfaces, i.e., in the par tial wetting regimes of 
the two individual surfaces [53|.15^. [ssj. or if one or both 
surfaces are covered by a thick wetting film [s^l- Such 
bridge formation may be relevant for colloid aggregation 
or flocculation of the particles [H, [57|] (for a summary of 
the corresponding experimental and theoretical work on 
these phenomena see Refs. [Uls^l)- For the sphere - pla- 
nar wall geometry relevant for the present experimental 
situation, theoretical studies (s^l predict that the bridg- 
ing transition can occur in the presence of thin wetting 
layers coating both surfaces. It is a first-order phase tran- 
sition and ends at a critical point. (Actually, these bridg- 
ing transitions are only quasi-phase transitions, because 
they involve, strictly speaking only a zero-dimensional 
volume [U, Is^l)- For a fixed distance between the wall 
and the sphere and fixed chemical potential, the posi- 
tion of this critical point is determined by the relation 
^ ~ i?, where R is the radius of the sphere (see Fig. d]). 
For small sphere radii the bridge configuration is unsta- 
ble, even for very small sphere-plane separations. On 
the other hand, bridging transitions are not possible for 
large sphere-plane separations, even if the sphere radii 
are large. The fiuid-mediated solvation force between the 
surfaces is very weak in the absence of the bridge and it 
is attractive and long-ranged if the capillary bridge is 
present. Moreover, for R/^ small its strength is propor- 
tional to the sphere- wall separation [s^, \EM j contrary to 
the case of two flat substrates [13] or to the sphere-sphere 
geometry [s^ . 



3. Critical Casimir forces for noncritical compositions 

For temperatures closer to the critical temperature 
the solvation force acquires a universal contribution 
due to the critical fluctuation of the intervening fluid 
which turns into the critical Casimir force. For a one- 
component fluid near gas- liquid coexistence = I-Lo(T) 
and confined between parallel plates it has been shown 
[4^ that at temperatures near the critical temperature 
Tc a smah bulk-like field h ~ = fi - no{T) < 0, 
which favors the gas phase, leads to residual condensa- 
tion and consequently to a critical Casimir force which, 
at the same large wall separation, is much more attrac- 
tive than the one found exactly at the critical point. The 
same scenario is expected to apply to binary liquid mix- 
tures, i.e., the Casimir force is expected to be much more 



attractive for compositions slightly away from the critical 
composition on that side of the bulk phase diagram which 
corresponds to the bulk phase disfavored by the confin- 
ing walls. This has been studied in detail in Ref. 
by using the standard field-theoretic model within mean- 
field approximation. These studies of the parallel plate 
geometry have been extended to the case of two spheri- 
cal particles of radius i? at a finite distance L [20] . The 
numerical results for the effective pair potential, as well 
as the results obtained by using the knowledge of the 
force between parallel plates and then by applying the 
Derjaguin approximation, valid for L <C i?, show that 
at T = Tc the dependence of the Casimir force on the 
composition exhibits a pronounced maximum at a non- 
critical composition. One expects that such a shift of the 
force maximum to noncritical compositions results from 
the residual capillary bridging and that the direction of 
the shift relative to the critical composition depends on 
the boundary conditions. If the surfaces prefer the a 
phase rich in species A, by varying the temperature at 
fixed off-critical composition c^, one observes that for 
small deviations |ca — c^^l <C cj^, the position of the max- 
imum of the Casimir force as function of temperature is 
almost unchanged, while the absolute value of the max- 
imal force increases considerably by moving away from 
to compositions ca < c'^. The overall temperature 
variation is, however, similar to that at c^, provided one 
stays sufficiently close to the critical composition. For 
compositions ca slightly larger than the critical compo- 
sition, Ca > c\, the critical Casimir force as a function 
of temperature is expected to behave in a similar way as 
for Ca < Ca, but the amplitude of the force maximum 
should be much weaker and should decrease for increas- 
ing Ca- At compositions further away from c^, i.e., off the 
critical regime, due to the small bulk correlation length 
the Casimir force is vanishingly small unless the afore- 
mentioned bridging transition is reached by varying the 
temperature. 

The case of a sphere against a planar wall has not 
been studied theoretically. However, we expect a similar 
behavior of the effective forces as the one for two spheres. 



III. EXPERIMENT 

A. The method: Total Internal Reflection 
Microscopy 

1. Basic principles of TIRM 

Total internal refiection microscopy (TIRM) is a tech- 
nique which allows one to determine the potential $ of 
effective forces acting on a single colloidal particle sus- 
pended in a liquid close to a planar substrate, with a force 
resolution down to the order of femto-Newton. The po- 
tential <I> is obtained from the probability distribution to 
find the surface of the particle at height z above the sub- 
strate, which is determined by monitoring the Brownian 
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FIG. 5: Data acquisition system (see main text for details). 
The green laser light generated by the optical tweezers is de- 
flected by a double prism into the microscope objective (in 
the figure represented as a grey vertical cylinder above the 
spherical particle) and it provides an optical potential which 
confines the spherical colloid laterally. The blue light which 
is scattered by this particle out of the evanescent field of in- 
tensity Iev{z) is collected by the same microscope objective, 
focused and then optically directed into the photon counter 
via a combination of prisms and mirrors (schematically rep- 
resented in the upper part of the figure). 

motion of the particle in the direction perpendicular to 
the substrate. In TIRM measurements this is achieved by 
creating an evanescent light field at the substrate-liquid 
interface which penetrates into the liquid. The intensity 
of the evanescent field varies strongly with the distance 
from the substrate. A single colloidal particle scatters 
light if it is illuminated by such a field. From this scat- 
tered intensity it is possible to deduce the position of the 
particle in the evanescent field, i.e., to determine z and 
its time dependence [H, [s^ . 

The basic experimental setup is presented in Fig. [S] A 
p-polarized laser beam (A = 473nm, P = 2 mW) is di- 
rected from below onto the interface between the bottom 
of a silica glass cell (a cuvette with a chamber to accom- 
modate a fluid film of thickness 200 /im) and the liquid 
containing the colloidal particle. The illumination angle 
9i (formed with the substrate normal) is larger than the 
critical angle 9c of total internal reflection. Due to total 
internal reflection, an evanescent wave penetrates into 
the medium with lower refractive index, here the liquid, 
and its intensity Icv{z) decays exponentially as a function 
of the distance z from the glass-liquid interface: 

/cv(^) = /ov(0)e-«^ . (16) 

The decay consta nt C defines the penetration depth 
which is given by [59| 



where A is the wavelength of the illuminating laser beam 
in vacuum, and rigiass and n^q are the refractive indices of 
the glass and the liquid, respectively. In our experiment 
(see, c.f., Sec. IIII B|) the critical binary mixture (liquid) 
has niiq = 1.384 whereas the silica glass (substrate) has 
f^giass — 1.464 (> niiq), resulting in 9c — 11° ■ A colloidal 
particle with a refractive index ncoU > J^iiq (in our ex- 
periment the polystyrene colloids have Ucow — 1.59) at a 
distance z away from the surface scatters light from the 
evanescent field. Within the well established data eval- 
uation model for TIRM intensity, the light scattered by 
the particle has an intensity Isc which is proportional to 
Isc oc Iev(,z) [5^ and therefore depends on the distance 
z. Care has to be taken in choosing parameters for the 
penetration depth and the polarization of the illuminat- 
ing laser beam in order to avoid optical distortions due 
to multiple reflections between the particle and the sub- 
strate, which would spoil the linear relation between /gc 
and Icv{z). In this respect, safe parameter regions are 
known to be small penetration depths C < 250nm and 
p-polarized illumination as used in the present experi- 
ment [60I . [6l| . As a result of this relation, the scattered 
light intensity Isc exhibits an exponential dependence on 
the particle-wall distance with exactly the same decay 
constant (C^^ = 200 ± 2nm in our experiment) as 
the evanescent field intensity lev'- 

/sc(z) = /o e-^^ (18) 

where the scattered intensity /q at contact z — Q depends 
on the laser intensity, the combination of refractive in- 
dices, and the penetration depth. As will be discussed 
below, the knowledge of Iq is important to determine the 
particle-substrate distance from the scattered intensity 
Isc- In principle in could be measured by the so-called 
sticking method [59|] according to which the particle is 
stuck on the substrate due to the addition of salt to 
the liquid in such a way as to suppress the electrostatic 
stabilization which normally repels the particle from the 
substrate. However, in the system we are interested in 
this is not practicable given the large concentration of 
salt (> 6mM) required to force the particle to stick to 
the surface and the compact design of the experimen- 
tal cell which limits the access to the sample. Instead, 
as described further below, we circumvent this problem 
by using a hydrodynamic method [i^l for the absolute 
determination of the particle-substrate distance. 

In Sec. |TT] we mentioned that, upon approaching the 
critical point of the binary liquid mixture, critical ad- 
sorption profiles form near the surfaces of the substrate 
and of the colloid. These concentration profiles induce a 
spatial variation of the refractive index, which deviates 
from the assumed steplike variation underlying Eqs. (jl6p 
and p?]) (see, e.g., Ref. [63]). Deviations from the func- 
tional form given by Eq. (|16p are pronounced if the corre- 
lation length ^ becomes comparable with the wavelength 
A of the laser light, which is not the case for the exper- 
imental data obtained here, for which A — 473 nm and 
^ < 100 nm (see, c.f., Figs.[Il[l5l imd^. 
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In a typical TIRM measurement run, the vertically 
scattered intensity Isc{t) = Isc{z{t)) (photons/s) is 
recorded by a photomultiplier connected to a single pho- 
ton counter [see Fig. [5] which counts the total number of 
scattered photons [see Fig. [DJa)] 



t+At 



(19) 



detected within a time interval At — 1 ms ^64] . The value 
of nsc{t) is then acquired with a frequency /samp = 250 Hz 
for a total duration tgamp — 15 min. The resulting set of 
data is then analyzed as described below in Sec. IIII A 21 
Consecutive intensity data Isc(t), i.e., nsc{t) acquired 
with a larger frequency /samp turn out to be strongly 
correlated in time. Accordingly, their acquisition does 
not contribute to the reduction of the statistical errors 
affecting the final estimate for the potential, as will be 
discussed in, c.f., Sec. lIVB"2l This observation motivates 
our choice /samp = 250 Hz. 

In addition to the detection optics, an optical tweezer 
is implemented in the TIRM setup [65] in order to be 
able to control the lateral position of the particle. The 
tweezer is created by a laser beam (Atweczcr = 532 nm) 
incident on the particle from the direction perpendicular 
to the substrate and focused by the microscope objective 
used also for the detection (see Fig. [5]). Via this tweezer 
it is possible to conveniently position the probe particle 
within the measuring cell and to restrict its lateral diffu- 
sion to a few microns so that the particle does not diffuse 
out of the field of view of the detection system. In addi- 
tion, the tweezer also exerts a light pressure |66j onto the 
particle, increasing significantly its effective weight (in 
the specific case considered here from ca. l.lkBTc/ i-iin to 
ca. T/cBTc/^m, see, c.f.. Sec. IIII A 31 and Refs. [67, ^ 
for details). In our experiment the tweezer is typically 
operated at a low power of P ^ 2 mW, but even at the 
largest power (Pmax — 25 mW) we used to trap and move 
the particle no effects of local heating, such as the onset 
of phase separation in the liquid, were observed due to 
the laser of the tweezer. 



2. Data analysis 

In order to determine the potential $ of the effective 
forces acting on the colloidal particle, one constructs 
a histogram out of the values of nsc{t) [see Eq. ([T^] 
recorded in the time interval tsamp, in such a way as to 
determine the probability distribution function Psc(^sc) 
for the particle to scatter rise = IscAt photons in a time 
interval At. Within the sampling time tsamp there are 
^samp/samp ~ N registration of counts. If N* is the num- 
ber of registrations which yield a certain count n*^, the 
probability p{nl^) of n*^ to occur is N* /N. By using 
Psc("-sc)dnsc = Pz{z)dz and Eq. (flSl) . this probability dis- 
tribution psc(''^sc) can be transformed into the probability 



for the particle-substrate distance z. In turn, in thermal 
equilibrium at temperature T, the probability Pz{z) is 
related to the particle-wall interaction potential ^{z) by 
the Boltzmann factor 



p,(z) = Cexp[-<I>(z)/(fcBT)], 



(21) 



where fesT is the thermal energy and C a normalization 
constant. Equation (|2ip holds because, due to the high 
dilution of the colloidal suspension, the single colloidal 
particle under observation does not interact with other 
particles. As a result, from the knowledge of Psc('T-sc) it is 
possible to determine $(z) up to an irrelevant constant 
related to C and to the overall normalization of Psc(''^sc)- 
For each bin of the histogram of Psc(''^sc), the correspond- 
ing distance z(nsc) is calculated via inversion of Eq. (11^ 



z(nsc) = -C ^ ln[nsc/(/oAt)] = 2:oxp(n-sc) - zq (22) 

where Zcxp("-sc) = — C^^lii'^sc is given in terms of ex- 
perimentally accessible quantities (i.e., tt-sc and C). This 
provides the position of the particle up to the constant 
Zq = — ln(/oAO the experimentally yet unknown 
position of the wall [Isc{z = 0) = /q]. In order to 
determine zq for all data sets, we have employed the so- 
called hydrodynamic method [62,] , which is based on the 
fact that due to hydrodynamic interactions the diffusion 
coefficient D of a colloidal particle at a distance z from 
the wall strongly depends on z. Moreover, near a wall 
the diffusion coefficient becomes also spatially anisotropic 
with the relevant value for TIRM measurements being 
D+, which refers to the diffusion occurring in the direc- 
tion perpendicular to the wall. Its spatial dependence 
can be expressed as 



D± = D^f{z/R) 



(23) 



where Doo = k-gT / [QTirjR) is the bulk diffusion coeffi- 
cient of a spherical particle of radius i? in a homogeneous 
fluid with viscosity 77 at temperature T. (For the water- 
lutidine mixture we use in our experiments, the value of 
rj has been measured in Ref. [70| as a function of tem- 
perature T and composition, with 77 ~ 2.09 x 10~'^Pas at 
T = 31°C and at the critical composition.) The reduced 
mobility function f{v) was calculated in Ref. [7l| and can 
be well approximated by 62]: 



f{v) 



6^2 + 2v 
6u2 + Qv + 2 



(24) 



Pz{z) = C'^sc(2)Psc("-sc(^)) 



(20) 



A plot of this theoretically predicted distance-dependent 
diffusion coefficient D± is shown in Fig. [T] A well es- 
tablished method [13] to determine the absolute particle- 
wall distance is to calculate the apparent diffusion co- 
efficient Z^app which is the weighted average of Da_{z) 
over the distances sampled by the colloidal particle, 
i.e., Dapp = /g°°dzDi(z)e-2f'^pz(z)//g°°dze-2'^'=p^(z), 
where the exponential factors in the numerator and de- 
nominator refiect the spatial dependence of /sc(^) [see 
Eq. and below]. This apparent diffusion coefficient 
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(a) (b) (c) 

FIG. 6: Data analysis: (a) Raw data for the total number of scattered photons nsc(t) detected at a certain time t within a 
time interval At = 1 ms from the single photon counter as a function of time, taken for a sampling time tsamp — 15 min. (b) 
Histogram calculated from this time series yielding the distribution function Psc(?^sc) of the number Usc of scattered photons, 
(c) The knowledge of the relation [see Eq. psp ] between the scattered intensity I^c — Tigc/At and the position z of the colloid 
in the evanescent field allows one to determine the probability distribution Pz{z) = C^sc{z)psc{nsc{z)) for the particle-substrate 
distance z from Psc(?^sc) and therewith by inversion of the Boltzmann factor the interaction potential. Further details are given 
in the main text. 
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FIG. 7: Diffusion coefficient D± of a spherical particle of ra- 
dius R moving perpendicular to a wall at a distance z [63,[z3l- 
Doo is the bulk diffusion coefficient. Note that for the present 
experimental conditions z/R < 1 and thus this distance de- 
pendence is pronounced. 



can be experimentally determined from the initial slope 
of the autocorrelation function C{dt) = (nsc(i)"-sc(i + 
6t))t of the scattering intensity fgll: 



-^app 



1 C'jO) 

e c{o) 



(25) 



where the prime ' denotes the derivative with respect to 
St. In order to determine zq one calculates the apparent 
diffusion coefficient i^app,caic on the basis of D±{z) and 
of the experimentally determined probability distribution 
Pz which is given by the parametric plot of C^scPsc('t-sc) 
as a function of Zoxp('T-sc) — zq upon varying Ugc [Eqs. (j20p 
and (|22p ]. with zq as the yet unknown position of the wall, 
and according to which the colloidal particle samples dis- 
tances. In turn, the value zq = zq can be determined by 



requiring that 



-Dapp.caic- A detailed description of 



B of Ref. ^6^) is primarily determined by the uncertain- 
tics of the particle radius [see Eqs. ((23)) and ([24]) ] and 
of the penetration depth [see Eq. (P5|) ]. Considering 
the experimental parameters and errors of our measure- 
ments, the resulting uncertainty in the particle-substrate 
distance z can be estimated to be ±30 nm for all plots 
shown in the following. 



this procedure can be found in Ref. [62| . The uncertainty 
in the determination of zq via this method (see Appendix 



We emphasize that it is sufficient to determine zq at a 
certain temperature in order to fix it for all the measured 
potential curves at different temperatures. Indeed the in- 
tensity /ev(0) of the evanescent field at the glass-liquid 
surface [as well as C~^, see Eg. pT)) ] depends on tempera- 
ture via the temperature dependence of the optical prop- 
erties of the glass and the liquid. In turn, this would im- 
ply a variation of the critical angle 9c with T, which was 
actually not observed within the range of temperatures 
investigated here. The intensity Iq oc /cv(0), which deter- 
mines Zq and which is recorded by the photomultiplier, is 
in principle affected by the temperature-dependent back- 
ground light scattering due to the critical fluctuations 
within the mixture (critical opalescence) . For the typical 
intensities involved in our experiment and for the temper- 
ature range studied, the contribution of this background 
scattering is actually negligible and, as a result, zq does 
not change significantly with temperature. The hydro- 
dynamic method, however, requires the knowledge of the 
viscosity rj of the mixture, which depends on tempera- 
ture and sharply increases upon approaching the critical 
point Wd\ due to critical fluctuations. These fluctuations 
might in addition modify the expression of D±{z). In 
order to reduce this influence of critical fluctuations we 
have chosen Tc — 3 K as the reference temperature for de- 
termining Zq, corresponding to a temperature at which 
no critical Casimir forces could actually be detected in 
the interaction potential. 
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3. Interaction potentials 

Under the influence of gravity, buoyancy, and tlic ra- 
diation field of the optical tweezer as external forces, the 
total potential $ of the colloidal particle floating in the 
binary liquid mixture, as determined via TIRM, is the 
sum of four contributions: 

$(z) = $0(2) + Goff2 + ^c{z) + ^offset . (26) 

In this expression $0 is the potential due to the electro- 
static interaction between the colloid and the wall and 
due to dispersion forces acting on the colloid; it is typ- 
ically characterized by a short-ranged repulsion and a 
long-ranged attraction. The combined action of gravity, 
buoyancy, and light pressure from the optical tweezer is 
responsible for the linear term GcffZ in Eq. (|^ (see, e.g., 
Ref. 66]). <&c(2) is the critical Casimir potential aris- 
ing from the critical fluctuations in the binary mixture. 
The last term 'l* offset is an undetermined, spatially con- 
stant offset different for each measured potential which 
accounts for the potentially different normalization con- 
stants of the distribution functions pi and Pz- While 
the flrst two contributions are expected to depend mildly 
on the temperature T of the fluid, the third one should 
bear a clear signature of the approach to the critical 
point. These expectations are supported by the experi- 
mental findings reported in Sec. IIVI The typical values 
of Gcff for the measurements presented in Sec. IIVI are 
Geff ~ 7.2fcBT/^m and GeS ~ lO.OfcsT/^m for the col- 
loids with diameters 2R = 2Aiim and 3.68/im, respec- 
tively [131 ■ Far enough from the surface, $0 and $c are 
negligible compared to the linear term and therefore the 
typical potential ^{z) is characterized by a linear increase 
for z large enough. Accordingly, upon comparing po- 
tentials determined experimentally at different tempera- 
tures, the corresponding additive constants ^offset, which 
are left undetermined by the TIRM method, can be fixed 
consistently such that the linearly increasing parts of the 
various coincide. However, it may happen that at some 
temperatures the total potential $ develops such a deep 
potential minimum that the colloid cannot escape from it 
and therefore the gravitational tail is not sampled. If this 
occurs the shift of this potential by a constant cannot be 
fixed by comparison with the potentials measured at dif- 
ferent temperatures. In order to highlight the interesting 
contributions to the potentials, the term GcffZ, common 
to all of them, is subtracted within each series of mea- 
surements and for all boundary conditions. Accordingly, 
the remaining part of the potential — displayed in the 
figures below — decays to zero at large distances. 

On the other hand, closer to the substrate, the (non- 
retarded) van-der-Waals forces contribute to ^o{z) with 
a term (see, e.g., Ref. Tab. S.5.b) 



ior $o,vdw(-z < R) 



-{A/6){R/z) to <i>o.ydw{z > 



R) 



-{2A/9){R/zy. The dependence of <^o, 



A 

<f>0,vdw(2) = - — 

6 



1 



1 



- ln(l + 2/6) 



(27) 



where A is the Hamaker constant and S = z/R. As 
6 increases, this term crosses over from the behav- 



R) on z is the same as the one of $0(2 <C ^,R) 
(see Eq. and therefore their relative magnitude 

is controlled by |$o,vdw(-2: < R)/^c{z < ^,R)\ = 
[^/(fcBT)]l/[6e(0)] ~ 2.4/(66(0)) for a critical point at 
T ~ 300 K and a typical value of the Hamaker constant 
A ~ 10-20 J. Taking into account that e(+,+)(0) ~ 2.5 
and 8(+__)(0) ~ 15 it is clear that in this range of dis- 
tances the critical Casimir potential typically dominates 
the (non-retarded) van-der-Waals interaction. Note, 
however, that for small values of z both of them be- 
come negligible compared to the electrostatic repulsion. 
For larger values of i? -C 2: ^ the Casimir force 
still dominates the dispersion forces, as discussed in de- 
tail in Refs. [l^, [5;^]. However, for the present ex- 
perimental conditions, the distance z is comparable to 
the bulk correlation length ^ and actually most of the 
data refer to the case £, ^ z < R, with the val- 
ues of the scaling variables ranging between z/^ ~ 
10 and z/R ~ 0.6 for distances at which the cor- 
responding potential is still experimentally detectable. 
Accordingly, with the above estimate for the poten- 
tial ratio, one can conservatively estimate |$ovdw(2; ~ 
0.6i?)/$c(2 ~ 1001 = [^/(fcBT)](0.2/|e(10j|) ~ 0.2 
where e(+,_)(10) ~ -e(+,+)(10) ~ 2.5. Retardation 
causes the van-der-Waals potential to decay as a function 
of z more rapidly than predicted by Eq. (|27p . This ad- 
ditionally reduces the contribution of ^o^vdw compared 
with (f>(7. In the analysis of the experimental data in 
Sec. IIVI we shall reconsider the Hamaker constant A for 
the specific system we are dealing with. 

In order to achieve the accuracy of the temperature 
control needed for our measurements we have designed a 
cell-holder rendering a temperature stability of ±5mK. 
This has been accomplished by using a fiow thermo- 
stat coupled to the cell holder with a temperature of 
(30.50 ± 0.01)°C functioning as a heat sink and shield 
against temperature fluctuations of the environment. In 
order to fine tunc the temperature the cell was placed 
on a transparent ITO (indium-tin-oxide) coated glass 
plate for a homogeneous heating of the sample from be- 
low. The voltage applied to the ITO coating was con- 
trolled via an Eurotherm proportional-integral-derivative 
(PID) controller for approaching the demixing tempera- 
ture. The controller feedback provides a temperature sta- 
bility of ±2mK at the position of the PtlOO-sensor used 
for temperature measurements. However, since the probe 
particle is displaced from the sensor by a few millime- 
ter, some additional fluctuations have to be considered. 
From the reproducibility of the potential measurements 
and from the relative temperature fluctuations of two in- 
dependent Ft 100 sensors placed on either side of the cell 
we inferred a ±5 mK stability of the temperature at the 
actual position of the measurement. The highly tempera- 
ture sensitive measurements were affected neither by the 
illuminating nor by the tweezing lasers due to moderate 
laser powers and due to low absorption by the probe par- 
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tide and by the surface. Although aU measurements were 
carried out upon approaching the demixing temperature, 
the Hght which is increasingly scattered in the bulk back- 
ground by the correlated fluctuations of the binary liquid 
mixture exposed to the evanescent field turned out to be 
negligible compared to the light scattered directly by the 
colloid. The effects of the onset of critical opalescence are 
significantly reduced by the fact that the illuminating op- 
tical field rapidly vanishes upon increasing the distance 
from the substrate and due to the still relatively small 
values of the correlation length. 



B. The binary liquid mixture and boundary 
conditions 

For providing the critical fiuctuations we have chosen 
the binary liquid mixture of water and 2,6-lutidine near 
its demixing phase transition. The bulk phase diagram 
of such a binary liquid mixture prepared at room tem- 
perature, ambient pressure, and sealed in a cell [73i] (con- 
stant volume) is reported in Fig. [51[a). It is characterized 
by a one-phase region (disordered phase) in which the 
two components form a mixed solution and which sur- 
rounds the closed loop of the two-phase region (ordered 
phase) in which these components segregate into a water- 
rich and a lutidine-rich phase. The first-order transition 
line delimiting the two-phase coexistence region, within 
which the two ordered phases form an interface, ends in a 
lower critical demixing point (LCP, see Fig. [Ha)) at the 
lutidine mass fraction Cj ~ 0.28 and the critical temper- 
ature Tc ~ 307.15 K [56|. The upper critical demixing 
point (UCP, see Fig. [5{a)) is located at high tempera- 
tures and therefore it occurs within the liquid phase only 
at pressure values above ambient ones. 

The choice of this specific binary liquid mixture as the 
critical medium is motivated by the fact that its prop- 
erties (bulk phase diagram, refractive index, viscosity, 
etc.) are known rather well and are documented in the 
literature, as this mixture has been extensively employed 
in the past for the study of phase separation per se or 
as a solvent of colloidal dispersions. A clear experimen- 
tal advantage is provided by the fact that the water- 
lutidine mixture at ambient pressure has a lower criti- 
cal point slightly above room temperature which can be 
conveniently accessed from the one-phase region by heat- 
ing the sample. Alternative binary liquid mixtures with 
a lower demixing point close to room temperature are 
formed by water and, e.g., 2,5-lutidine (Tc ~ 13°), 2,4- 
lutidine {Tc ~ 23°C), triethylamine (T^ ~ 18°C), and 
n-butoxyethanol {Tc ~ 49°C). The addition of a suitable 
amount of a third component to some of these binary liq- 
uid mixtures allows one to shift the critical temperature 
of demixing basically at will. For example, by adding 
3-methylpyridine to the presently used mixture of 2,6- 
lutine and water it is possible to gradually shrink the 
coexistence loop in Fig. [HJa), causing an upwards shift 
of Tc of as much as 30°C. Analogously, even though 3- 



methylpyridine is always miscible with normal water, it 
exhibits a coexistence loop if mixed with heavy water, 
with a lower critical point at Tc ~ 39° C. The addition 
of normal water to this demixing binary liquid mixture 
of heavy water and 3-methylpyridine causes an upwards 
shift of Tc until the coexistence loop disappears at a dou- 
ble critical point with T^'^p = 78°C The variety 
of available substances and the tunability of the critical 
temperature by suitable chemical additions allow one to 
generate critical Casimir forces for a conveniently wide 
range of temperatures. 

In our experiment the mixture is prepared under nor- 
mal conditions (room temperature, ambient pressure) 
and then it is introduced into the sample cell which is 
afterwards sealed with Teflon plugs in order to hinder 
the mixture from evaporating. Although we have no con- 
trol on the resulting pressure p of the mixture, the fact 
that the cell is not hermetically sealed and that small air 
bubbles might be trapped within it should keep p very 
close to its ambient value. Within the limited range of 
temperatures we shall explore in the experiment, possible 
pressure variations are not expected to lead to substan- 
tial modifications of the phase diagram (e.g., shifts of the 
critical point) compared to the ones at constant pressure 
or volume. The order parameter (j) for the demixing phase 
transition can be taken to be the difference between the 
local concentration Ci,(x) {mass fraction) of lutidine in 
the mixture and its spatially averaged value cl- Accord- 
ingly, a surface which preferentially adsorbs lutidine is 
referred to as realizing the (-I-) boundary condition for 
the order parameter given that it favors (j) > 0, whereas 
a surface which preferentially adsorbs water leads to the 
(— ) boundary condition. 

The experimental cell containing the binary liquid mix- 
ture and the colloid is made up of silica glass. Depend- 
ing on the chemical treatment of its internal surface, one 
can change the adsorption properties of the substrate so 
that it exhibits a clear preference for either one of the 
two components of the binary mixture. In particular, 
treating the surface with NaOH leads to preferential ad- 
sorption of water (— ), whereas a treatment with hexam- 
ethyldisilazane (HMDS) favors the adsorption of lutidine 
(+) [3 1 as we have experimentally verified by compar- 
ing the resulting contact angles for water and lutidine on 
these substrates. As colloids we used polystyrene parti- 
cles of nominal diameter 2R — 3.69 /im and 2.4/im, the 
latter possessing a rather high nominal surface charge 
density of 10/iC/cm'^. Size polydispersity of these par- 
ticles are 2% and 3%, respectively, corresponding to ca. 
±70 nm. (These nominal values are provided in the data- 
sheets of the company producing the batch of particles, 
see Ref. [3^ for details.) The adsorption properties of 
polystyrene particles in a water-lutidine mixture have 
been investigated in Refs. [t^, [t^, with the result that 
highly charged (> 3.8/iC/cm^) colloids preferentially ad- 
sorb water (highly polar) whereas lutidine is preferred 
at lower surface charges. Even though we did not in- 
dependently determine these adsorption properties, the 
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(a) (b) 

FIG. 8; Bulk phase diagram of the binary liquid mixture of water and 2,6-lutidine (dimethylpyridine C7H9N) [s^ . [t^ I at 
constant volume. The relevant thermodynamic variables are the temperature T and the mass fraction cl of lutidine in the 
mixture. Open symbols in panel (a) refer to actual experimental data [t^. The schematic side view of a vertical sealed cell 
filled with the binary liquid mixture is shown by the insets (i) and (ii) of panel (a) for thermodynamic states outside and 
inside of the coexistence loop, respectively. In (ii) W and L indicate the water- and the lutidine-rich phase, respectively. The 
mixture separates into the lutidine-rich and the water-rich phase within the two-phase coexistence area encircled by the solid 
first-order transition line. At the lowest and highest points (LCP and UCP, respectively) of this line the demixing transition is 
continuous. The detailed phase diagram in the vicinity of the lower critical point is shown in panel (b), together with the typical 
thermodynamic paths (dashed vertical lines) experimentally investigated here within the gray region as well as for cl = 0.2. 



(particle, substr.) 


colloid diam. 27?: 


substrate treat.: 


3.69/im 


2.4^m 


HMDS 


(+,+) 


(-,+) 


NaOH 


(+,-) 


(-,-) 



TABLE I: Experimental realization of all possible (particle, 
substrate) symmetry-breaking boundary conditions, where 
(-I-) indicates the preferential adsorption of lutidine and ( — ) 
the preferential adsorption of water. The treatments of the 
substrate affect only its surface properties. 



the surface properties of the colloidal particles we used 
in the experiment. In order to rule out a possible degra- 
dation of the colloid during the experiment, we verified 
the reproducibility of the observed effects after each data 
acquisition. 



IV. RESULTS 



results of Refs. [tsI . [76j and the corresponding nominal 
values of the surface charges of the colloids employed in 
our experiment suggest that the polystyrene particles of 
diameter 2i? = 3.69 /J,m [2R = 2.4 have a clear pref- 
erence for lutidine (-I-) [water )]. A posteriori, these 
presumed preferential adsorptions are consistent with the 
resulting sign of the critical Casimir force observed ex- 
perimentally. Depending on the surface treatment of the 
cell and the choice of the colloid one can realize easily all 
possible combinations of (particle, substrate) boundary 
conditions (see Tab.|T]). 

For a given choice of the particle-substrate combina- 
tion with its boundary conditions and for a given con- 
centration Cl of the mixture we have determined the in- 
teraction potential $ (see Eq. (PS|) ) between the colloid 
and the substrate as described in the previous subsection, 
starting from a temperature T below the critical point in 
the one-phase region and then increasing it towards that 
of the demixing phase transition line at this value of c^. 
It might happen that, as a result of leaching, the water- 
lutidine mixture slowly (i.e., within several days) alters 



In, c.f.. Figs. [TU1 - [T51 and [T71 we report the experimen- 
tally obtained interaction potentials $ as functions of the 
particle- wall distance z for various values of the tempera- 
ture T, both at critical (Figs. [TUHT5|) and off-critical con- 
centrations (Fig. fTT)) . In all the cases presented, the grav- 
itational and offset parts of the potentials [see Eq. (p6|) ]. 
which turn out to be de facto independent of the tem- 
perature T , have been subtracted in such a way that the 
resulting potentials vanish for large values of z. How- 
ever, those potentials, for which the gravitational tail 
could not be sampled (see, e.g., Figs . [T^ and [T7|) . cannot 
be normalized like the others by this requirement. 

Depending on the concentration of the mixture, two 
qualitatively different behaviors are observed, which are 
discussed in Sec. IIV Bl for cl — cl and in Sec. IIV Cl for 
Cl ^ c'j^. However, in the next subsection we first discuss 
the experimental results for the potentials measured far 
away from the transition line and the comparison of them 
with theoretical predictions. This provides important in- 
sight into the effective background forces to which the 
critical Casimir forces add upon approaching the critical 
point. 



16 



A. Non-critical potentials 

In all the cases reported here, sufficiently far from the 
transition line one observes a potential which appears 
to consist only of the electrostatic repulsion between the 
colloid and the substrate and which can be fitted well by 



$o(z) = fceTe-'^^^-^-) 



(28) 



where is the Debye screening length and Zos the value 
of the distance z at which $o(-z = Zcs) — k^T. (zcs is ex- 
pected to depend, inter alia, on the surface charge and on 
the radius of the colloid.) For the potential in, c.f., Fig.fTOl 
which corresponds to Tc — T = 300 mK, a fit of k yields 
12nm ± 3nm, which is compatible with the esti- 



K 

mate n' 



10 nm derived from the standard expression 



K = \/e2X;iPi/(eiiq(0)fcBT) (see, e.g., Ref. [7g|), where 
e is the elementary charge, eiiq(O) the static permittivity 
of the mixture (see below) , and pi the number density of 
ions assumed to be monovalent and estimated in Ref. [TSj 
for the dissociation of a salt-free water-lutidine mixture. 
Within the range of distances z sampled in our experi- 
ment there is no indication of the presence of an attrac- 
tive tail in $0: which on the other hand is generically 
expected to occur due to dispersion forces, described by 
a potential as given in Eq. (j27p . In order to compare 
this experimental evidence with theoretical predictions, 
below we shall discuss the determination of the Hamaker 
constant A in Eq. (|?7|) on the basis of the dielectric prop- 
erties of the materials involved in the experiment. The 
relation between them is provided by (see, e.g., Ref. [T^ I 



A{z)^%tY^' ^ 

n— ^ 



^liq ^coW ^liq 



Rn{z) 

(29) 



where the permittivities e of the various materials are 
evaluated at the imaginary frequencies zwn, with LOn — 
2TTkBTn/h = 71 X 2.5 x IQi^'rad/s at T ~ BOOK. (Note 
that the imaginary part of the complex permittivity e{ijj) 
as a function of the complex frequency u vanishes on the 
imaginary axis Rew = (zll-) The factor Rn{z) accounts 
for retardation and, neglecting the fact that in the three 
different media light propagates with different velocities 
(i.e., for Rn assuming ecoU — euq — £coii) it takes the 
form Rn = (1 -I- r„)e~''". The ratio r„ = 2r„(z)/a;~^ 
quantifies the relevance of retardation: heuristically, a 
thermally fluctuating electric dipole within, e.g., the 
glass generates an electric field which travels at least 
a distance z across the liquid, taking a minimal time 
T,i{z) = z/{c/ ■y/eiiq(iti;„)), before inducing an electric 
dipole within the colloid. Such an induced dipole, in 
turn, generates an electric field which travels back to the 
original dipole and interacts with it. However, such an 
interaction is reduced by the fact that the original dipole 
has a lifetime uj~^ and might have decayed during the 
minimal time 2r„(2;) it takes the electric field to do the 



roundtrip which is the case for r„ = 2t„(z)/u;~^ > 1. 
The prime in Eq. indicates that the contribution of 
the static permittivities n = is to be multiplied by 1/2 
(see, e.g., Ref. [72]), resulting in the term An=o- Within a 
first approximation, in Eq. (j29[) the retardation factor i?„ 
does not affect those terms with ?■„ < 1, corresponding to 
Rn — 1 for them, while it exponentially suppresses those 
terms with r„ ^ 1. Accordingly, only the former con- 
tribute significantly to ^„>o and retardation is accounted 
for by summing in A„>o only the terms corresponding to 
^ (c/ y^e{iu!n))/{2z), whcrc s{iuJn) ^ niiq (see, c.f., 
dash-dotted line in Fig. [5] for n = 1, . . . , 4) is the refrac- 
tive index of the liquid. The expression in Eq. (|29[) is 
valid generically for a film and only in the non-retarded 
regime i?„ ~ 1 for the sphere-plate geometry. However, 
even for the latter geometry an estimate of the order of 
magnitude of the effects of retardation can be inferred 
simply by restricting the sum in Eq. ([29| to the values of 
n corresponding to rn ^ I, so that i?„ ~ 1. 

The parameters of our experiment, i.e., z > 0.1 /im 
and niiq ~ 1.38 yield t^„ < 1.1 x 10^^ Hz, i.e., n < 4 at 
T ~ 300 K, with ujn=i,...,4 in the infrared (IR) spectrum. 
The contribution An=o of the zero-frequency mode is ac- 
tually subject to screening by the salt in the liquid so- 
lution, characterized by the Debye screening length 
which also controls the exponential decay of the elec- 
trostatic contribution to $o- This means that An=o 
is not a constant but acquires a z dependence. This 
is accounted for by a multiplicative correction factor 
Ro{z) = (1 -I- 2Kz)e~^'^^ multiplying An=o as given by 
Eq. (|^ : we note that i?o(0) ^ I > Ro{z). A parametric 
representation of the permittivity e^^{iL^J) of pure water 
can be found in Table L2.1 in Ref. [72]. In the part 
of the spectrum over which the sum in Eq. (|29p runs, 
for polystyrene the permittivity ecoii(*t^) is actually al- 
most constant £coU = ^^coii — 2-53 (see, e.g.. Table L2.3 
in Ref. [t^ for a parameterization). For lutidine and 
silica, the parameters which characterize the correspond- 
ing permittivities ei^{iuj) and sgis^ssii'-^) are summarized in 
Table 1 of Ref. [tJ- In order to calculate the dielectric 
permittivity euqiiuj) of the homogeneous water-lutidine 
mixture on the basis of ew, £l, and the lutidine volume 
fraction — 0.25, one can use the Clausius-Mossotti 
relation, as explained in Ref. [t^- The resulting permit- 
tivities are reported in Fig. [HI With these elements at 
hand and within the approximation discussed above one 
can calculate an upper bound to the value of A„>o which 
accounts for retardation, finding 

An>oiz > 0.1 fim) < 0.06 /cb^ ~ 0.025 x lO-^^J. (30) 



In addition, from Eq. 
eters one finds 

Ar^ 



and from the material param- 



0.46 fceT 



(31) 



resulting in a screened contribution A^"q (z) = 
An^oRo^z), with ~ 12nm, which is actually neg- 
ligible compared to the electrostatic potential $0 [see 
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FIG. 9: Permittivities e{iuj) of the materials relevant for our 
experiment as functions of the frequency u; [rad/s] on a loga- 
rithmic scale. As on the left, from bottom to top, we report 
the curves corresponding to pure lutidine, polystyrene (EcoIi), 
silica (sgiass), a water- lutidine mixture with a volume concen- 
tration <I)L — 0.25 of lutidine (euq), and pure water. These 
curves are based on the parameterizations and on the mate- 
rial properties reported in Refs. [t^, [z3- The four vertical 
lines indicate the frequencies {uJn}n=i,...4 which within our 
approximation enter into the determination of the Hamaker 
constant in Eq. (|29|l . (Note that lme{Lj) — if, as in the 
present case, Retj = [73].) 



Eq. 4"o(2)/$o(z) < A«=oe-""- < 6x 10-4 (d., 

Tab. [IT] for typical values of k and Zcs)- Taking into ac- 
count Eq. ((27l) . the contribution of A^^"q to $o,vdw(^) 
becomes comparable to the electrostatic one only for 
z/R < 10"'*, which is well below the range z/R > 0.1 
actually investigated in our experiment. On the other 
hand, for z ~ 0.1 /im, ^^^"q — 0.004 /cb^T and due to its 
exponential decay with z, we expect this contribution to 
be negligible compared to the value Anyo{z ^ 0.1 /im) in 
Eq. pO|) . In turn, this value is considerably smaller than 
the typical one A ~ 10~^° J we have used in Sec. IIII A 31 
in order to compare dispersion forces with the critical 
Casimir potential. Accordingly, the conclusion drawn 
there that the latter typically dominates the former is 
confirmed and reinforced by the estimate for A given 
here. 

The values just determined for An=Q and A„>o are 
meant to be estimates of their orders of magnitude given 
that a detailed calculation which properly accounts for 
retardation (following, e.g., Ref. [73]) and for possible in- 
homogeneities in the media, especially within a binary 
liquid mixture, goes beyond the present scope of a quali- 
tative comparison of theoretical predictions for the back- 
ground forces with the actual experimental data. In this 
respect, a detailed determination of the permittivities of 
the specific materials used in our experiment would be 
crucial for an actual quantitative comparison of this con- 
tribution to $ with the experimental data. With all these 
limitations, the theoretical calculation discussed above 
yields A{z > 0.1 /im) ~ A„>o(z > 0.1 /im) < 0.06 fceT. 



If one insists on fitting the experimental data for the 
background potential ^o{z) by including the contribu- 
tion of the dispersion forces as given by Eq. (|27p in addi- 
tion to a possible overall shift A$offsct, one finds values 
for the Hamaker constant A which vary as function of 
the range of values of z which the considered data set 
refers to. This might be due to the fact that the statisti- 
cal error affecting the data increases at larger distances 
or due to an incomplete subtraction of the gravitational 
contribution, which might bias the result. In particular, 
in the range < z < 0.3 fim we focus on data for the 
potentials which have been measured experimentally for 
the largest temperature deviation from the critical point 
and which are smaller than 6 fcsT. The choice of this 
latter value results from a compromise between avoid- 
ing the increasing statistical uncertainty due to the poor 
sampling of the sharply increasing potential and having 
a sufficiently large number of data points left at short 
distances, where ^o{z) is not negligible. The resulting 
parameter values for the four experimentally measured 
potentials are reported in Tab. |TT1 The resulting values 
of A are compatible with a rather small Hamaker con- 
stant, in qualitative agreement with the previous theo- 
retical analysis. The combined estimate of the screen- 
ing length is somewhat larger than anticipated from the 
analysis of one of the potentials [see after Eq. (|28p ] and 
results in = (14 ± 4) nm, again in agreement with 
independently available experimental data [TSj . In or- 
der to highlight the presence of dispersion forces in this 
system, here masked by the strong electrostatic repul- 
sion, one would have to increase the salt concentration of 
the solvent in order to reduce significantly the screening 
length which then provides access to smaller particle- 
substrate distances. However, we emphasize that a de- 
tailed and quantitative study of these background forces 
is not necessary in order to identify the contribution of 
critical Casimir forces to the total potential and it is 
therefore beyond the scope of the present investigation. 



B. Critical composition 

1. Experimental results 

For the binary liquid mixture at the critical composi- 
tion we have estimated (after data acquisition) the criti- 
cal temperature Tc as the temperature at which anoma- 
lies in the background light scattered by the mixture 
in the absence of the colloid and due to critical opales- 
cence are observed and visual inspection of the sample 
displays an incipient phase separation. The value deter- 
mined this way has to be understood as an estimate of 
the actual value of the critical temperature of the water- 
lutidine mixture and it is used for the calibration of the 
temperature scale, which is shifted in order to set Tc to 
the nominal value Tc = 307.15 K reported in the litera- 
ture (see, e.g., Ref. ^). Note, however, that depend- 
ing on the different level of purity of the mixture, pub- 
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BC 


Fig. 


Tc - T[K] 




A^offBctfteT] 


K ^ [nm] 


ZcB [nm] 


(-,-) 


|lo| 


0.30 


0.2 ±0.1 


0.14 ±0.08 


10.5 ±0.5 


113 ± 1 


(+,-) 


|11| 


0.90 


0.2 ±0.1 


0.2 ±0.1 


17 ± 1 


90 ±3 


(+,+) 


|12| 


0.20 


0.05 ±0.03 


0.06 ±0.04 


15.9 ±0.5 


85 ± 1 


(-,+) 


|13| 


0.31 


0.0 ±0.2 


0.0 ±0.1 


13 ± 1 


153 ±2 



TABLE II: Fit parameters for the non-critical potentials 4>non-cr for four boundary conditions and with the gravitational part 
subtracted, <I>non-cr — "l'o,vdw(2) ± A$offsct ± fcBTe~'°'^~^'==', where "l>o,vdw(2) is given by Eq. ([27ll . The values reported here 
correspond to 95% confidence intervals for the parameters. 
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FIG. 10: Effective interaction potential ^{z) between a wall 
and a spherical particle of radius 1.2 /xm immersed in a water- 
lutidine mixture at the critical concentration c^, as a func- 
tion of the distance z from the wall and for various values of 
the temperature in the one-phase region (T < Tc) [l^]- The 
gravitational and the offset contribution to the potential [see 
Eq. (|26p ] have been subtracted. The set of solid and dashed 
lines, which are barely distinguishable on this scale, corre- 
spond to the theoretical predictions (see the main text for de- 
tails). The potentials reported here refer to the (— , — ) bound- 
ary conditions (other cases are reported in Figs. [TT] and I12p 
and show that an increasingly attractive force contributes to 
the total potential upon approaching the critical tempera- 
ture, i.e., upon decreasing AT = Tc — T. Here Tc is the 
nominal value of the critical temperature, corresponding to 
the anomaly in the background scattering, which signals the 
onset of critical opalescence in the sample. Only the data to 
the right of the vertical dotted line are considered for the com- 
parison with the theoretical predictions (see the main text). 

lished experimental values of Tc are spread over the rang e 
306.54 307.26 K (see, e.g., the summary in Ref. [83|). 
Due to the difficulties in determining the absolute value 
of the critical temperature, with our experimental setup 
only temperature differences are reliably determined and 
the actual critical temperature of the mixture might dif- 
fer slightly from the nominal value Tc- We shall account 
for this fact in our comparison with the theoretical pre- 
dictions. 

Close to Tc critical opalescence is expected to occur. It 



FIG. 11: Interaction potential ^(z) as in Fig . 1101 for the 
(±, — ) boundary conditions and R ~ 1.85 /im [37|. An in- 
creasingly repulsive force contributes to the total potential 
upon approaching the critical temperature. The set of solid 
and dashed lines, which are barely distinguishable on this 
scale, correspond to the theoretical predictions (see the main 
text for details). Only the data to the right of the vertical 
dotted line are considered for the comparison with the theo- 
retical predictions (see the main text). 

is indeed ultimately observed upon heating the mixture 
towards the critical temperature, leading to an increase 
in the background light scattering due to the correlated 
fluctuations in the mixture. Even though this might in- 
terfere with the determination of the interaction poten- 
tial $ via TIRM, within the range of temperatures we 
have explored at the critical concentration, the enhance- 
ment in the background scattering is actually negligible 
compared to the light scattered by the particle. 

In Fig. [TO] we present the interaction potentials <!> as 
a function of the distance z for that choice of colloidal 
particle and surface treatment which realizes the (— , — ) 
boundary condition (see Tab.|T|. As discussed above, for 
AT = Tc - T = 0.30K, the potential consists only of the 
electrostatic repulsion (see Eq. (^5)) ). Upon approach- 
ing the critical point an increasingly deep potential well 
gradually develops, indicating that an increasingly strong 
attractive force is acting on the particle. At the smallest 
AT we have investigated, i.e., AT = 0.12K, the resulting 
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FIG. 12: Interaction potential ^(z) as in Figs. llOland llll for 
the (+, +) boundary conditions Wq. As for ( — , — ) boundary 
conditions (see Fig. llOf) . upon approaching the critical tem- 
perature an increasingly attractive force contributes to the 
total potential. 



potential well is so deep that the particle hardly escapes 
from it. In view of the small temperature variation of 
ca. 180 mK, the change of ca. IO/cb^ in the resulting po- 
tential is remarkable. This very sensitive dependence on 
T is a clear indication that in the present case critical 
Casimir forces are at work. In the case of Fig. [10] the 
maximum attractive force acting on the particle is about 
600 fN. 

According to the theoretical predictions, one expects 
the critical Casimir force to be repulsive for asymmet- 
ric boundary conditions (+, — ) or (—,+). In our ex- 
periment we can easily realize the — ) BC by chang- 
ing the colloidal particle surface from preferentially ad- 
sorbing water (— ) to preferentially adsorbing lutidine 
without any additional surface treatment of the 
cell (see Tab. HI. The interaction potentials <& for this 
case are reported in Fig. [TT] As for the (— , — ) bound- 
ary condition, sufficiently far from the critical point (i.e., 
AT = Tc-T = 0.90 K), the potential consists only of the 
electrostatic repulsion contribution [see Eq. (pS)) ]. How- 
ever, upon approaching the critical point the repulsive 
part of the potential curves shifts towards larger values 
of the distance z, indicating that an additional repulsive 
force is acting on the colloid. It is possible to make this 
force attractive again by treating the surface of the cell so 
that its preferential adsorption changes from water (— ) 
to lutidine (-I-), so that then (+,+) boundary condition 
is realized. As Fig. [T^] shows, the resulting potentials 
show indeed an attractive part the qualitative features of 
which resemble those of the case with the (— , — ) bound- 
ary condition, reported in Fig. (TO] Note that the depth 
of the potential in Fig. [T^] corresponding to AT = 0.05 
is so large that the gravitational part (which has been 
subtracted) cannot be sampled by the particle and there- 



FIG. 13: Interaction potential <l>(z) as in Figs.[Tni[ni andfT^ 

for the (—,-!-) boundary conditions. As in the case of the 
(-I-, — ) boundary conditions (see Fig. Illf) an increasingly re- 
pulsive force contributes to the total potential upon approach- 
ing the critical temperature. The set of solid and dashed 
lines, which are barely distinguishable on this scale, corre- 
spond to the theoretical predictions (see the main text for 
details). Only the data to the right of the vertical dotted 
line are considered for the comparison with the theoretical 
predictions (see the main text). 



fore the position of this potential curve along the vertical 
axis cannot be fixed. If, with the same (-I-) surface of the 
cell, one changes again the colloidal particle from prefer- 
entially adsorbing lutidine (-I-) to preferentially adsorb- 
ing water (— ), we can experimentally realize the (— ,+) 
boundary condition (see Tab. |T| for which a repulsive 
critical Casimir force is expected. The resulting poten- 
tial is reported in Fig. I13l and shows the same qualitative 
features as the one in Fig. 1111 with an increasingly repul- 
sive force which builds up upon approaching the critical 
point. 



2. Comparison with theory 

The experimental data reported in the previous sub- 
section can be compared with the theoretical predictions 
presented in Sec. Ill A[ which are expected to be valid 
for S = z/R <C 1 (Derjaguin approximation). In the 
experimental setting corresponding to Figs. [TU] and 1131 
R — 1.2/^m whereas z < 0.3/im and z < 0.8/im, re- 
spectively, so that S < 0.25 and S < 0.67. In Figs. [11] 
and [T^ one has R ~ 1.85/im with z < 0.8 and z < 0.3, 
respectively, corresponding to S < 0.43 and 5 < 0.16. 
Accordingly the Derjaguin approximation is expected to 
provide a sufficiently accurate description of the experi- 
mental data, possibly apart from those at larger values of 
z in Figs. [TT] and I13[ the corresponding potential values 
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of which are anyhow neghgibly small. In order to ex- 
tract from the measured potential only the part which 
is due to the critical Casimir force we focus on that 
range of distances z for which the electrostatic contri- 
bution ^q{z) (see Eqs. (|26l) and ((28|) ) as measured far 
from the critical point (i.e., for AT = 0.30K in Fig. [TOl 
and AT = 0.90K in Fig. [TT|l is actually negligible, using 
as a criterion |$o(^)| ^ O-SfceT. This latter choice also 
avoids additional complications due to possible changes 
induced by critical fluctuations in the electrostatic con- 
tribution $0(2^) upon approaching the critical point. Ac- 
cordingly, for a quantitative comparison with the theo- 
retical predictions we consider only data corresponding 
to z > 0.12/zm in Fig. [TO] [(-,-) BC], z > O.l^m in 
Fig. [m [{+, -) BC], and z > 0.16/.tm in Fig. [TO] [(-, +) 
BC], excluding in each case the data set corresponding 
to the largest value of AT, which has been used to define 
^o{z). Unfortunately, the number of data points which 
satisfy this condition in the case of Fig. [T^] is quite lim- 
ited for providing a basis for a reliable analysis; therefore 
we do not process these corresponding data. 

The strength of the critical Casimir force depends 
strongly on the deviation Tc — T from the critical point 
via the bulk correlation length ^ ~ |Tc — T\~^ . Accord- 
ingly, even a small systematic error in the experimental 
determination of T and can result in sizeable discrep- 
ancies in the comparison between measured potentials 
and theoretical predictions. Statistical variations and a 
possible drift of the temperature during the acquisition 
of the data, which are kept within 5 mK by the temper- 
ature controller used in our experimental apparatus [37| , 
are similarly important. In order to circumvent parts 
of these problems we compare the experimental data, se- 
lected by the aforemention criterion (see the vertical lines 
in Figs. [TOl El and [TO)) , for a certain boundary condition 
and for the six temperatures Ti = Tc~ ^Ti, i — 1,. . . , 6, 
closest to Tc with the theoretical prediction for ^ci^; ^) 
provided by Eqs. (O and ([TO]). For each Ti the val- 
ues of the correlation length and of a possible residual 
offset ^offset,! are determined in such a way as to opti- 
mize the agreement between $c(^; Ci) + 'i'offset.i and the 
corresponding experimental data set. A drift of the tem- 
perature during the acquisition of the data might affect 
the value of the correlation length resulting from this 
procedure. As we shall see below, even if present, this 
possible drift does not strongly affect the final estimate 
for the correlation length amplitude ^0 in Eq- ©, the 
uncertainty of which is dominated by the systematic un- 
certainty of the theoretical predictions stemming from 
finite-size extrapolations of the Monte Carlo data. 

The data set (Ti,^i) is then fitted with the theoreti- 
cally expected power-law behavior given by Eq. ([2|) . This 
is carried out by fixing the universal critical exponent 
ly ~ 0.630 to its best known theoretical value while de- 
termining the non-universal amplitude ^q*^*"* and the value 



^t^^ [A] 


Ref. (year) 


Method 


2.0 ±0.2 


[83] (1972) 


static LS 


2.92 ±0.19 


[83] (1972) 


dynamic LS, linewidth 


2.7 ±0.2 


[84] (1987) 


static LS 


2.3 


[85] (1993) 


specific heat 


2.1 


[86] (1997) 


critical adsorption 


1.98 ±0.04 


[82] (2006) 


specific heat 



TABLE III: Experimental estimates of the non-universal cor- 
relation length amplitude ^0 for the water-lutidine mixture 
at the critical concentration. In light scattering (LS) ex- 
periments the bulk correlation length ^{T) is determined by 
measuring the wave-vector (static) or frequency (linewidth, 
dynamics) dependence of the scattered intensity. A fit of 
£,{T) to the expected algebraic behavior [Eq. ^] yields the 
value ^0- Alternatively, ^0 can be obtained on the basis of 
the measured value of the non-universal amplitude A'^ which 
characterizes the divergence of the specific heat at constant 
pressure Cp(r — > 0^) ~ (^"*'/a)r"" and the theoretical [s^] 
or experimental value of the universal amplitude ratio 
Q"*" = yl"'"5op/(MfcB), where M is the molar mass of the mix- 
ture and p its mass density at the critical point [HI . [l^ . A 
careful theoretical analysis of experimental data for the criti- 
cal adsorption profiles also leads to an estimate for ^0 [86l ]. 
Comparing the experimentally measured potentials to the 
theoretical predictions for the critical Casimir contribution 
we obtain the estimate reported in, c.f., Eq. (|34|) . 

ri^*'' of the critical temperature from the data set (T^, ^j): 

Here we assume that the temperatures Ti are sufficiently 
close to Tc so that ^ is described correctly by its leading 
power-law behavior. The resulting value of ^q^*'' can then 
be compared with the available independent experimen- 
tal estimates reported in Tab. lIIIi providing a check of the 
consistency of the experimental data with the theoreti- 
cal predictions. In spite of the scattering of the available 
experimental data, which might be due to different con- 
ditions of the mixture (such as contaminations or slightly 
different concentrations) or to different systematic uncer- 
tainties of the various approaches, all the estimates are 
within the range 

^^"''P^ = 2.3 ±0.4 A, (33) 

estimated via a least-square fit of the data [s^, [H, [13, 
Issl . Ispl in Tab. IIIII (The experimental value quoted in 
Ref. |37l| refers to the estimate of Ref. [s^.) The limited 
set of temperatures which have been investigated exper- 
imentally (apart from one far away from Tc and used for 
fixing the background potential, six different values for 
each set of boundary conditions) does not allow us to de- 
termine simultaneously and reliably the exponent from 
the experimental data. 

In comparing the experimental data with the theoret- 
ical predictions we have to take into account the un- 
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certainties which affect both of them. As discussed in 
Sec. Ill Al the currently available theoretical predictions 
within the Derjaguin approximation are affected by a 
20% systematic uncertainty, clearly visible in Fig. [21 for 
the amplitude of the scaling function of the Casimir po- 
tential $c- 

As far as the experimental data are concerned, the sys- 
tematic uncertainties — which are the ones most relevant 
for the comparison — concern (i) the particle-wall dis- 
tance z, which can be determined by the hydrodynamic 
method up to Az = ±30nm (see SecIIISand Ref. ^) 
and (ii) the absolute temperature scale T and, in par- 
ticular, the value of the critical temperature Tc- As de- 
scribed above, in order to cope with the uncertainty in T 
we opted for an indirect determination of the associated 
correlation length from the best fit of the experimental 
data with theoretical predictions. 

In addition to these systematic uncertainties, there are 
statistical errors associated with the way the potential <I> 
is determined via TIRM. The number N{nsc) of counts 
during the sampling time tsamp, registered in each bin 
of size Arise, centered around n^c and forming the inten- 
sity histogram pscl'T.sc) = A^("sc)/(A'^tot AtIsc) reported in 
Fig. [6fb), is subject to statistical fluctuations AN{nsc) 
which affect the estimate of Psc("-sc)i Pz{z), and there- 
fore of the potential A'^tot — /sampisamp is the total 
number of counts in the time series of Uscit) of dura- 
tion tsamp, from which the histogram of psc(fisc) has been 
constructed. One expects that these statistical fluctua- 
tions are relatively more important for those bins which 
are less populated, i.e., for smaller N{nsc)- In terms of 
the distance z of the colloid from the wall, they corre- 
spond to values which are less frequently sampled during 
the Brownian motion of the particle under the influence 
of the potential i.e., to larger values of the poten- 
tial. This can be seen directly from the potentials re- 
ported in Figs, m [m m [H and, c.f., [HI in which the 
experimental data are more scattered very close to the 
wall and far from it, whereas the sampling of the po- 
tential $ is particularly accurate around its minimum. 
In order to evaluate the statistical uncertainty associ- 
ated with each data point of the potential, ideally one 
should construct the histogram of Pscirisc) based on sev- 
eral different realizations of the time series nsc{t) and 
then analyze the statistical properties of this ensemble 
of plots. Alternatively, one might evaluate the autocor- 
relation time tcorr of nsc(t) [e.g., from a detailed study 
of the autocorrelation function C{6t), see the text before 
Eq. (HH)]. Assuming that the number A^in(?T-sc) of sta- 
tistically independent counts in a bin of the histogram 
is given by Mn(nsc) = A^('^sc)/(/samptcorr) and assum- 
ing that the statistics of the counts is Poissonian, the 
associated relative statistical fluctuation is relate d to 
the number of counts by AA'in(nsc) = -^in(^sc) and 
induces a statistical uncertainty A^{z{nsc))/{kBT) — 

ANinirisc) / Nin{nsc) = ^in^^^("-sc) for the value of the 
potential $ at the position zijigc) corresponding to the 
scattered number of photons rigc- However, in the com- 



parison between the experimental data and the theoret- 
ical predictions, the statistical error in the former is ex- 
pected to be negligible compared to the systematic un- 
certainty in the latter and therefore we do not proceed 
to a detailed evaluation of the statistical error associated 
with the data points reported in Figs. [TUl El [HI [131 , 
and, c.f.,[T71 Actually, a good estimate of the magnitude 
of the statistical error can be inferred from the scatter of 
the experimental data points relative to a smooth curve 
interpolating each potential. 

In order to discard those data which are affected by 
large statistical fluctuations, we consider for the compar- 
ison with the theoretical prediction only data fulfilling 
^(z) < 6kBT and z < z,nax where Zmax = 0.3 /xm for 

Fig.iini 
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for Fig. [T3l For the (— ,— ) boundary conditions, in 
Fig. [TUl we report the comparison between the theoret- 
ical prediction and those experimental data which have 
been selected as explained above. The solid and dashed 
lines, barely distinguishable, correspond to the predic- 
tions given by Eq. (fT2l) by using for 0(_._) — ©(+.+) 
the scaling functions described by the solid and dashed 
lines, respectively, in Fig. [21 (For each temperature ATi 
the experimental potentials have been shifted vertically 
by the amount — ^offost.i determined previously as the 
best fit parameter.) The corresponding values £,i of the 
correlation length are reported in Fig. [141 together with 
the resulting best fit based on Eq. which leads 

to the least-square estimates ^q^'"* = 1.7 ± O.lA and 
Ari*^*' =Tc- Ti^*^ = 52 ± 10 mK for the sohd curve and 
= 1.7±0.1A and Ari^*^ = 65±7 mK for the dashed 
one. Taking into account the systematic uncertainty of 
the scaling functions for the Casimir potential, we arrive 

at the combined estimate ATi^*'' = 60± 15 mK, i.e., the 
value of Tc determined experimentally is actually higher 
than the value Tc^^^ resulting from the comparison with 
the theoretical predictions. In addition, this comparison 
allows one to estimate the correlation length ^, for which 
no independent experimental estimate is presently avail- 
able. According to Fig. [TH one has 20nm < ^ < 40nm, 
so that for the range 0.12/im < z < 0.3/xm of distances 
this translates into the ranges 3<x = 2;/^<15 and 

6 < u = (1 - T/Tf '^/'^(z/C^*^*^/'' ^ 70 of the scal- 
ing variables x and u (see Fig. [5]). In order to be able 
to test prominent features of the theoretically predicted 
scaling function such as the occurrence of a minimum for 
Umin — 0.5, one has to reach ^ > 180 nm, i.e., one must 
get still closer to the critical point (AT < 6 mK) than it 
was possible in the present experiment. 
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FIG. 14: Correlation length ^ as a function of the temper- 
ature deviation AT — Tc — T from the experimentally lo- 
cated critical temperature Tc for the water-lutidine mixture 
at the critical concentration. The sets of data points are ob- 
tained by optimizing the agreement between the experimental 
data from Fig. [TD] and the corresponding theoretical predic- 
tions for the Casimir potential for the ( — , — ) boundary con- 
ditions (Fig. [21 see the main text for details). The upper 
(lower) set of points has been determined by using as the 
scaling function = 0(_|.__|.) of the Casimir potential 

the one reported by the dashed (solid) line in Fig. [21 The 
dashed and solid curves are the best fit of the corresponding 
data sets based on Eq. p2p with the corresponding values 
Ari^*' = Tc — Ti*''' indicated by the vertical dashed and 
solid lines, respectively, i.e., ^ diverges at Ti*"'-* corresponding 

to Ar(r = Ti*^'') = Tc - Ti*"> = Ari"*'. 



For the (+, — ) boundary conditions, in Fig. [TT] we 
report the comparison between the theoretical predic- 
tion and those experimental data which have been se- 
lected as explained above. As in Fig. (TUl the solid 
and dashed lines correspond to the predictions based 
on Eq. by using for 0(+^_) = ©(-,+) the scaling 
functions indicated as solid or dashed line, respectively, 
in Fig. 121 (As in Fig. [TOl the experimental potentials 
have been shifted vertically for each temperature ATi 
by the amount — $offcst,i determined previously as the 
best fit parameter.) The corresponding ensuing values 

of the correlation length are reported in Fig. [151 to- 
gether with the resulting best fit based on Eq. ([32l) , which 
leads to the least-square estimates ^q*^*-* = 1.8±0.lA and 
^jXfit) ^ _ y^(fit) ^ 234 ± 2mK for the solid curve 
and = 1.9 ± O.lA and At]*^*^ = 235 ± 2mK for 

the dashed one. The final combined estimate of ATi^*-*, 
which takes into account the systematic uncertainty of 
the amplitude of the theoretical prediction of 0(+^_), is 

therefore ATi^*-* = 235 ± 3 niK. The correlation lengths 
reported in Fig. [TSl are in the range 20nm ^ £. ^ 95nm. 
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FIG. 15: Same as Fig. 1141 obtained by optimizing the agree- 
ment between the experimental data from Fig. [11] and the 
corresponding theoretical predictions for the critical Casimir 
potential for the (-I-, — ) boundary conditions (Fig. [2}. 



The corresponding ranges of distances z depend on the 
temperature AT (see Fig. [111) so that the experimental 
data cover the scaling variable ranges 3.4 ^ x = z/^ < 17 

and 7 < u = (1 - T/Tf *')V"(-z/fr')'^" ^ 85 (see 
Fig. [2]). 

For the (— ,+) boundary conditions, in Fig. [T3l we 
report the comparison between the theoretical predic- 
tion and those experimental data which have been se- 
lected as explained above. As in Figs. [TOl and (TTJ the 
sets of solid and dashed lines, barely distinguishable, 
correspond to the predictions based on Eq. (fT2|) by us- 
ing for = the scaling functions indicated 
as solid or dashed fine, respectively, in Fig. [2] (As 
in Figs. \TU\ and [TTl the experimental potentials have 
been shifted vertically for each temperature AT^ by the 
amount — $offcst,i determined previously as the best fit 
parameter.) The corresponding values S,i of the correla- 
tion length, which can be inferred from this comparison, 
are reported in Fig. [111 together with the resulting best 
fit based on Eq. ([32l) . which leads to the least-square es- 
timates C^^*^ = 2.55 ± 0.25A and At]*^*' = Tc - t]*^*' = 
14 ± 11 mK for the solid curve and ^^^^^ = 2.7 ± 0.2A 
and Ari***^ = Tc - T^:^^^ = 10 ± 9 mK for the dashed 
one. The final combined estimate of ATi'^*^ which takes 
into account the systematic uncertainty of the ampli- 
tude of the theoretical prediction of 6(_^+), is therefore 

(fit'] 

ATc — 14 ± 11 mK. The correlation lengths reported 
in Fig. [ini are in the range 28 nm ^ C ^ 66 nm with the 
corresponding ranges of the distances z depending on the 
temperature AT (see Fig. [T5|) so that the experimental 
data cover the scaling variable ranges 7<x — z/^<67 

and 22 < li = (1 - T/Ti^'V^^'C-^/CrV^'' ^ 790 (see 
Fig. [2]). 
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comparison with (-,+) BC 
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FIG. 16: Same as Figs. [14] and 1151 obtained by optimizing 
the agreement between the experimental data from Fig. [13] 
and the corresponding theoretical predictions for the critical 
Casimir potential for the ( — , +) boundary conditions (Fig.[2]|. 



The experimental data reported in Fig. [T3] (which 
Fig.dHrefers to) have actually been acquired by an exper- 
imental setup which makes use of an improved tempera- 
ture control compared to the one used during the acquisi- 
tion of the data reported in Figs. [TOl dH (H and[I7| This 
upgrade of the setup is characterized by a better temper- 
ature stability and allows one to determine Tc with higher 
accuracy [st} . 

As far as the value of is concerned, taking into ac- 
count the values reported above for the (— , — ) and {+, — ) 
boundary conditions, one obtains the combined estimate 
C^^*^ = 1.8 ±0.2 A [3^ which is in very good agreement 

with the experimental value Co'"^''^ reported in the first 
line of Tab. IIIII (and quoted in Ref. [33|)- It is interesting 
to note that the principal source of error in these esti- 
mates of ^0 is actually the systematic uncertainty in the 
theoretical predictions, which turns out to be more sig- 
nificant than the statistical or possible systematic exper- 
imental errors, such as the one due to possible variations 
or fluctuations of the temperature occurring during the 
measurement. Actually, due to the pronounced depen- 
dence of the theoretical predictions on the temperature 
via the correlation length, these variations should result 
in averaged effective values of the correlation length, most 
probably affecting the overall amplitude, i.e., the value 
of The estimate of ^^^'^ = 2.6 ± 0.3 A based on 

the data for the (— , boundary conditions agrees with 
the estimate reported in the third line of Tab. IIIII but it 
is larger and not quite compatible with the former one. 
This might be due to larger systematic errors in the lat- 
ter or due to possibly different conditions of the mixture 
employed in the experiment (e.g., purity or possible con- 
tamination by leaching). The combined estimate which 



accounts for the results of our analysis is therefore 

4^'^ = 2.2 ±0.6 A, (34) 

which is compatible and similar to the experimental value 
^(cxp) j.gpQj.|;g(j gq (j33p })Yit carrying a larger, mainly 
systematic, uncertainty. This agreement is particularly 
significant if one takes into account the fact that Eq. (jM)) 
combines results obtained from different experimental 
conditions (different particles and different surface treat- 
ments) , interpreted on the basis of the available theoret- 
ical predictions. It is worthwhile to point out that one 
can equally reverse the line of reasoning given above. One 
can adopt the point of view that the correlation length is 
a bulk property which is determined by independent bulk 
measurement yielding Eqs. ^ and This way ^ is 

not a fit parameter but an input which fixes the theo- 
retical prediction for the critical Casimir potential com- 
pletely. Since ^l^^^^^ equals Ca^^\ this implies that with 
this fixed input the theoretical predictions yield the solid 
and dashed lines in Figs.[Tni[Ill and [131 In this sense one 
can state that the remarkable agreement between theory 
and experiment observed in the figures occurs without 
adjusting parameters. 

C. Noncritical composition 

In the sense of renormalization-group theory the devi- 
ation of the bulk concentration cl from its critical com- 
position c'l represent the second relevant scaling field be- 
sides the reduced temperature. In the language of the 
Ising universality class the field conjugate to this de- 
viation plays the role of a bulk magnetic field. (Here 
we do not discuss that actually two linear combinations 
of the conjugate field and of the reduced temperature 
form the appropriate orthogonal scaling fields.) In this 
sense, as already mentioned in Sec. IIIBI measurements 
along thermodynamic paths of varying temperature at 
fixed off-critical compositions probe the dependence of 
the critical Casimir forces on another, equally important 
scaling variable. Corresponding theoretical predictions 
have been derived for the parallel-plate geometry |88l| as 
well as for two adjacent spherical particles [l^, l20j|. one 
expects that a similar behavior holds for the present ge- 
ometry of a sphere near a planar surface. For a composi- 
tion of the binary liquid mixture far away from its critical 
value the critical Casimir forces become negligible. For 
small deviations from the critical lutidine mass fraction 
c1 and for the (+,+) BC (so that both surfaces pref- 
erentially adsorb lutidine) the temperature variation of 
the Casimir force upon approaching the two-phase coex- 
istence line near the critical point from the mixed phase 
depends on whether cl is larger or smaller than c|^. If 
is slightly less than the critical composition, cl < cl, the 
Casimir force is expected to behave similarly as along 
the critical composition, i.e., it should exhibit a mini- 
mum before reaching the temperature T^,{cl) at which 
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FIG. 17: Effective interaction potential $(2) between a wall 
and a particle of diameter 2R = 3.69 fim, immersed in a water- 
lutidine mixture at mass fraction cl = 0.2 < function 
of the distance z and for various values of the temperature 
close to the demixing phase boundary [37|. This system corre- 
sponds to the (+, +) boundary conditions. The gravitational 
and offset contributions to the potential [see Eq. (I26p ] have 
been subtracted. The abrupt formation of a narrow minimum 
of the potential close to the wall upon increasing the tempera- 
ture is interpreted as the formation of a liquid bridge between 
the particle and the wall (see the main text). 



phase separation occurs [T^{cl — c^) = Tc\. However, 
the depth of the critical Casimir potential is expected to 
be considerably deeper away from the critical composi- 
tion (see Fig. 8 in Ref. [131 )• For lutidine mass fractions 
slightly larger then the critical value, cl > cl, as func- 
tion of temperature the Casimir force is expected to vary 
similarly as in the case < but much weaker. 

These expectations are in agreement with the obser- 
vations made in our experiment for the (-I-, BC and 
for several values of lutidine mass fractions in the range 

0. 26 < Cl < 0.32 [not shown]. For cl close to the 
measured potentials between the wall and the colloidal 
particle look similar to those obtained for the critical 
composition. For cl much smaller than the potentials 
are similar to the ones shown in Fig.[T7](see, e.g.. Fig. 6.9 
in Ref. [8§|, corresponding to cl ~ 0.25). 

For the values of cl further away from the system is 
no longer near criticality. Therefore the critical Casimir 
force ceases to influence effective interactions between the 
colloidal particle and the wall. However, for the (-I-, 4-) 
BC and for temperatures allowing for phase separation, 

1. e., in the present case above the critical temperature, 
one expects a bridging transition to occur at composi- 
tions Cl < c^. If the lutidine concentration is below 
its critical value the conjugate bulk field favors the 
water-rich phase whereas for the (-1-, +) BC the confin- 
ing surfaces prefer the lutidine-rich phase. At a single 
wall this competition gives rise to wetting phenomena 



and in confined geometries to condensation phenomena. 
As discussed in detail in Sec. IHB[ if a bridge of the con- 
densing phase connects two adjacent spheres immersed in 
the binary liquid mixture, there is an attractive wetting- 
induced interaction that pulls the sphere together. We 
expect a similar scenario to occur for the present geome- 
try, i.e., if a wetting bridge of the phase which is favored 
by both surfaces is formed between the spherical particle 
and the planar wall. 



Indeed, the measurements carried out for lutidine mass 
fractions cl < 0.2 indicate that such a bridge formation 
takes place. In Fig.[T7]we plot the measured particle- wall 
interaction potentials at several temperatures near but 
below the temperature of demixing T^{cl = 0.2), which 
could not be located with sufficient accuracy. Moreover, 
the temperature scale in Fig. [17] has not been calibrated 
with a reference temperature, so that only temperature 
changes are significant. As one can see for temperatures 
between 307.31 and 307.36 K, the potentials are very well 
described by Eq. ((26)l which accounts for the electrostatic 
and gravitational contributions only. However, upon fur- 
ther increasing the temperature by 20 mK we observe a 
markedly different behavior of the potentials. Suddenly 
the interaction potentials are shifted towards smaller val- 
ues of distances z. Also the shape is changed in that the 
potentials exhibit a narrow and deep minimum. This 
sudden shift of the potential well towards the surface in- 
dicates the onset of an attractive interaction between the 
particle and the wall. A further slight increase in tem- 
perature gives rise to an even stronger shift of the po- 
tential minimum towards the wall. This phenomenon is 
observed only for lutidine mass fractions smaller than 
the critical value, i.e., on that side of the phase diagram 
where the mixture is poor in the component that is pref- 
erentially adsorbed by both surfaces. This behavior is 
in stark contrast to critical Casimir forces which vary 
gradually as function of the thermodynamic variables. 
On the other hand, for the (— ,— ) boundary conditions 
and Cl ^ cl the resulting potentials do not differ qual- 
itatively from those shown in Fig. [TOl But by further 
increasing the concentration of lutidine to cl — 0.4 one 
observes the sudden formation of a narrow and deep po- 
tential well upon increasing the temperature (see the po- 
tential corresponding to cl — 0.3, 0.32, and 0.4, reported 
in Fig. 6.11 of Ref. 89]). These observations are in agree- 
ment with the theoretical concepts described in Sec. Ill Bl 
The effective potentials associated with the formation of 
a bridge, formed by the phase coexisting with the bulk 
phase, between the particle and the substrate are the- 
oretically expected to exhibit hysteresis upon changing 
the temperature back towards its start value. However, 
in the present experiment, only rather weak hints for this 
hysteresis have been observed and actually no convincing 
evidence for it could be produced. 
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V. SUMMARY, CONCLUSIONS, 
PERSPECTIVES, AND APPLICATIONS 

A. Summary 

We have presented a detailed account of the exper- 
imental and theoretical investigations of the effective 
forces acting on spherical colloidal particles of radius R 
close to a substrate and immersed in a near-critical bi- 
nary liquid mixture, shortly reported in Ref. ^37*1 . Based 
on total internal reflection microscopy (Figs. O [HI and[7]) 
our main experimental flndings are the following: 

1. Upon raising the temperature T of the binary liquid 
mixture of water and lutidine at its critical concen- 
tration towards its lower critical point Tc of demix- 
ing (see Fig. [5]), an attractive or repulsive force act- 
ing on the colloidal particle arises gradually. 

2. This effective force is attractive if the surfaces of the 
colloid and of the substrate display preferential ad- 
sorption of the same component of the mixture (see 
Figs.[TUland ll2p . whereas it is repulsive in the cases 
of opposing preferences (see Figs. [Til and fT3)) . This 
contribution to the total effective force (compare 
Sec. IIV Al and Fig.lHl) is negligible at temperatures 
a few hundred mK away from Tc and it increases 
significantly upon approaching it. As experimen- 
tally verified, these so-called critical Casimir forces 
can be reversibly switched on and off by changing 
the temperature. 

3. If the concentration of this binary liquid mixture is 
close to but not equal to the critical one we have 
observed a behavior which is qualitatively similar 
to the one observed for the mixture at its critical 
concentration. In the close vicinity of the critical 
point there is no experimental evidence for the oc- 
currence of wetting phenomena. 

4. If the concentration of the binary liquid mixture 
differs significantly from the critical one and both 
surfaces exhibit the same preferential adsorption 
for that component of the mixture which is disfa- 
vored in the bulk, we observe the abrupt formation 
of a narrow and deep potential well (see Fig. [T7|) 
upon approaching the phase boundary of first-order 
demixing. 

The experimental observations 1, 2, and 3 can be con- 
sistently interpreted in terms of the occurrence of the 
critical Casimir effect in near-critical mixtures, whereas 
observation 4 can be understood in terms of the forma- 
tion of a bridgelike configuration of the segregated phases 
(see Fig. U) . For mixtures at the critical concentration 
it is possible to quantitatively compare the measured po- 
tentials with the corresponding theoretical predictions for 
the contribution of critical Casimir forces (see Eqs. (fTU)) 
and (|12p as well as Figs. [2] and [3]), derived within the 
Derjaguin approximation (see Fig. [T]) and for the range 



of distances within which electrostatic forces are negli- 
gible (see the solid fines in Figs. [H El and [131). The 
correlation length ^, as determined from the comparison 
between the experimental data and the theoretical pre- 
dictions, follows rather well the theoretically expected 
universal power-law behavior (see Figs. [HI [151 ^^^d \W\i 
and the associated non-universal amplitude is in agree- 
ment with previous independent experimental determina- 
tions for this specific binary mixture [compare Eq. (|34[) 
to Eq. dm and see Tab.lm]. 

The same critical Casimir forces, investigated here by 
using a water-lutidine mixture, are expected to act on a 
colloid immersed in any binary liquid mixture close to 
its demixing point (or in any fluid close to its gas-liquid 
critical point) and in the vicinity of a substrate. While 
the values of non- universal parameters, such as Tc and 
^0, depend on the specific mixture, the resulting criti- 
cal Casimir force is described by a material- independent, 
universal scaling behavior (see Eq. Q) and scaling func- 
tion (see Eq. ^ and Fig. [3l for small particle-substrate 
separation) which depends only on whether the adsorp- 
tion preferences of the particle and of the substrate are 
equal [(+, +), (-, -)] or opposite [(+, -), (-, +)]. 



B. Discussion 

The experimental observations summarized above 
might contribute to the understanding of the reversible 
aggregation of a dilute suspension of colloidal particles 
immersed in a water-lutidine mixture close to its demix- 
ing point, which has been the subject of several exper- 
imental studies since it was first observed in 1985 [s^ 
(for a review see Ref. ) ■ The formation of pre- wetting 
layers around the particles was first invoked as a possi- 
ble explanation for this phenomenon. Later on it was 
experimentally demonstrated that aggregation might ac- 
tually occur in a region of the phase diagram which ex- 
tends too far from the two-phase coexistence line and 
from the wetting transition to be possibly related to 
pre- wetting phenomena ^76j. (However, no aggrega- 
tion was observed for mixtures at the critical concentra- 
tion [t^.) Among the possible different mechanisms (see, 
e.g., Ref. fri\ for a summary) which might contribute to 
explain this flocculation, also critical Casimir forces have 
been invoked theoretically, as summarized and discussed 
in Refs. [l^, [l^,!!^. In particular, the experimental ob- 
servation (besides for 2,6-lutidine and normal water as 
solvent reported also for colloids dispersed in mixtures 
of 3-methylpyridine, heavy, and normal water [tsI . [79j 
or 2-butoxyethanol and normal water [sO]) that floccu- 
lation phenomena arc enhanced near but off the criti- 
cal point, at compositions which are slightly poorer in 
the component preferentially adsorbed by the colloids 
than the critical one, matches with the fact that also 
the critical Casimir forces attain their maximum values 
there. Although flocculation involves the interaction of 
many colloidal particles and therefore is a many-body 
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phenomenon (which can also be inter pre ted as a genuine 
phase transition in a ternary mixture [8l|), some features 
such as the experimentally determined asymmetry of the 
aggregation line with respect to the critical concentration 
(see, e.g., Refs. [H, [t^) can be qualitatively accounted 
for by the behavior of the effective interaction among two 
colloids, mediated by the near-critical solvent. This prob- 
lem was theoretically investigated in Refs. [I9i, 20. 39i| 
within various approximations. In this context, our ex- 
perimental study of the interaction between a single col- 
loidal particle and a substrate suggests that the actual 
magnitude of the forces due to the critical Casimir effect 
and due to the formation of a bridgelike configuration are 
large enough to play an important role also in aggrega- 
tion phenomena. This has been demonstrated recently 
on patterned substrates [iol [ooj. However, more quanti- 
tative corresponding statements require additional ded- 
icated studies of many-body effects. We note that, de- 
pending on the specific physical and chemical properties 
of the colloidal suspension under consideration, various 
mechanisms might be at play in determining its aggrega- 
tion, especially for charged colloids, for which screening 
effects or even field-induced pha se separation of the mix- 
ture might be predominant [77l [Qlj. In this respect, re- 
cent experimental studies of the structure factor of such 
an aggregating colloidal suspension via synchrotron small 
angle X-ray scattering (SAXS) 79] might provide impor- 
tant insight into the physical mechanisms at play in that 
phenomenon. 

At the critical composition, the present experimental 
study detected the occurrence of critical Casimir forces 
in a range of distances z which corresponds to a scal- 
ing variable x = z/^ ^ 1. In this limit, some of the 
qualitative features of the scaling function of the force 
(such as the occurrence of a maximum for (-I-, bound- 
ary conditions — see Fig. [3]) have not been probed. Ac- 
tually, in this limit the associated potential for (+, ±) 
boundary conditions is very well described by an expo- 
nential function $c(^)/(fcB7') ^ 2nA±{R/^) exp{-z/^) 
[see Eqs. and HH)]. A clear signature of the collec- 
tive nature of such an interaction is the fact that its range 
is set by the correlation length ^. On the other hand, the 
functional form of this dependence on z is quite generic 
and actually is common to interactions of rather differ- 
ent nature (e.g., electrostatic ones - see Eg. ipS)) - the 
range of which is set by the screening length k^^). For 
a relatively small correlation length ^ ~ the critical 
Casimir interaction and the electrostatic repulsion have 
the same range and, depending on the specific values of 
the parameters, one of them might dominate over the 
other. Especially in this case one expects an interplay 
between these two effects due to the fact that the order 
parameter profile develops inhomogeneities on a length 
scale ^ which is comparable to the typical length 
which characterizes the electrostatic screening in the ho- 
mogeneous medium. Taking into account that due to the 
segregation of the mixture the ions have different solu- 
bilities in water and lutidine, the enhancement of one 



of these two components close to the confining surfaces 
might result in a change of the screening of the electro- 
static interaction compared to the case of a non-critical 
homogeneous medium. Analogously, as dispersion forces 
depend sensitively on the spatial structure of the dielec- 
tric media forming the system and on the associated per- 
mittivities e(tj, x), the inhomogeneities which build up in 
the medium upon increasing ^ might affect significantly 
the background van-der-Waals contribution $o,vdw to 
the total potential compared to the estimate we gave for 
a homogeneous mixture [see Eqs. ^7} and (|29p ]. In con- 
trast to the critical Casimir effect, however, a quantita- 
tive analysis of the interesting interplay between critical 
fluctuations and dispersion/electrostatic forces necessar- 
ily requires the knowledge of several system-specific prop- 
erties such as the actual spatially varying composition of 
the mixture and the resulting permittivity e(a;, x). Some 
of these properties might be inferred experimentally via, 
e.g., surface plasmon spectroscopy of the binary mixture 
close to the substrate. The comparison of the experimen- 
tal evidences presented here with the theoretical predic- 
tions has not generated an actual need to account for the 
possible interplays mentioned above. Presumably they 
result into effects which are negligible in the range of 
variables explored in our experiment and within our ex- 
perimental accuracy. This might not be the case for dif- 
ferent choices of the particle and the mixture for which, 
e.g., dispersion forces and therefore their possible modifi- 
cation due to critical fluctuations might be more relevant 
than in the system investigated here. 

In order to compare the experimental data with the 
theoretical predictions we have inferred the bulk correla- 
tion length ^ from the experimental data for the critical 
Casimir potential (see Sec. IIVBI) . Reversing the line of 
argument, ^ can be inferred on the basis of the deviation 
AT from the critical temperature and of the knowledge 
of the non-universal amplitude which has been deter- 
mined by independent bulk experiments. However, one 
could also determine the scaling functions of the criti- 
cal Casimir potentials on the sole basis of experimental 
data, without need of any additional theoretical informa- 
tion. This requires the experimental determination of the 
actual correlation length ^ corresponding to each temper- 
ature, which can be accomplished, e.g., by studying the 
light scattered by the mixture sufficiently far from the 
colloid and the substrate, probing the behavior of bulk 
fluctuations. Such an independent experimental deter- 
mination of the scaling function would provide an addi- 
tional, valuable test of the theoretical predictions. 



C. Perspectives and applications 

Suitable chemical treatments of the surfaces in con- 
tact with the binary liquid mixture can be used in order 
to control the strength of their preferential adsorptions 
(i.e., the corresponding surface fields, see Sec. Ill A]) and 
therefore the resulting critical Casimir force. The exper- 
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imental data presented here have been consistently inter- 
preted in terms of the theoretical predictions correspond- 
ing to strong preferential adsorption. However, in view 
of possible applications of these effective forces, it is also 
important to study in detail both experimentally (ooj and 
theoretically [H, \^ 13] cases in which such a pref- 
erence is weaker or spatially modulated in a controllable 
fashion via suitable chemically patterned substrates. In 
the latter case and depending on the symmetry of the 
pattern, the resulting critical Casimir force acting on the 
colloidal particle acquires a lateral component in addi- 
tion to the normal one investigated here [40l . [90j . This 
lateral force, as the normal one, is characterized by a uni- 
versal scaling behavior and its range is again set by the 
correlation length ^, such that it can be switched on and 
off by controlling the distance from the critical point. In 
addition, the force turns out to be rather sensitive to de- 
tails of the imprinted chemical structure, e.g., the striped 
pattern considered in Refs. 4fl"2fl|- ^ proper theoretical 
analysis of the critical Casimir potential enables one to 
infer from the experimental data knowledge about such 
details even if they could not be determined by indepen- 
dent means, such as atomic force microscopy !40]. 

The lateral Casimir force might also find applications 
in colloid rheology. Consider, e.g., a dilute suspension of 
(+) colloids exposed to a suitably fabricated substrate 
which has an adsorption preference smoothly varying 
along one direction from (— ) to (+), such that it changes 
appreciably on the scale of the radius of the colloid. For 
sufficiently small values of the correlation length, the col- 
loids diffuse isotropically along the substrate. However, 
upon approaching the critical point, the lateral Casimir 
force associated with the spatial gradient of the prefer- 
ential adsorption (i.e., of the corresponding surface field) 
adds a deterministic drift to this diffusion process, re- 
sulting in a transport of colloids along the surface of the 
substrate. The direction of the flow will be reversed by 
changing the preferential adsorption of the colloid from 
(-I-) to (— ), which can be exploited as a reversible se- 
lection mechanism. In this context, the critical Casimir 
force acting on a micrometer-sized colloid exposed to a 
substrate with a modulated adsorption preference on the 
the scale of some hundred micrometers has been recently 
studied experimentally 87]. 

Topographical modulations of the surface of an other- 
wise chemically homogeneous substrate can also be used 
to control the direction of the total force acting on a 
similar substrate fv?\ or on a colloidal particle exposed 
to it. Additionally, chemical patterning or geometrical 
deviations from spherical symmetry of the colloidal par- 
ticle (e.g., ellipsoidal colloids) result in a critical Casimir 
torque [35, 36] if the particle is exposed to a substrate. 
Combining all these features one should be able, e.g., to 
control reversibly via minute temperature changes the 
orientation of such colloids exposed to a striped sub- 
strate. 

The critical Casimir force acting on a colloidal par- 
ticle close to a plate fluctuates in time due to the fact 



that it originates from time-dependent critical fluctua- 
tions. In the present analysis we focused on the mean 
value of such a force and on the associated averaged 
potential <I>c- However, as explained in Sec. pil A[) . 
TIRM naturally provides a measurement of the time- 
dependent sphere-plate distance z{t) which, in turn, can 
be used to determine the correlation time of the crit- 
ical Casimir force and how its expected algebraic tem- 
poral singularity builds up upon approaching the critical 
point. (Some aspects of this dynamical behavior are dis- 
cussed in Refs. [H, [qI, ^^.) This critical slowing down 
of the critical Casimir effect can in principle be exploited 
in order to control the resulting dynamics of the colloidal 
particle. 

In contrast to interactions which typically act among 
colloids (such as electrostatic and dispersion forces), the 
critical Casimir force is characterized by a pronounced 
temperature dependence. This fact can possibly be 
exploited in order to control via minute temperature 
changes the phase behavior and aggregation phenomena 
in colloidal dispersions in the bulk or close to those chem- 
ically structured solid surfaces which find applications in 
the fabrication of nano- and micrometer scale devices. 
Not only the range of the resulting interaction can be 
controlled but also its sign and spatial direction. This 
can be typically achieved by surface treatments and it 
does not require (as it does, e.g., for dispersion forces) 
substantial changes or a fine tuning of the bulk properties 
of the materials which constitute the immersed objects 
and the mixture itself. These properties could be used, 
e.g., in order to neutralize the attractive quantum me- 
chanical Casimir force responsible for the stiction which 
brings micro-electromechanical systems to a standstill. If 
these machines would work not in a vacuum but in a liq- 
uid mixture close to the critical point, the stiction could 
be prevented by tuning the critical Casimir force to be 
repulsive via a suitable coating of the various machine 
parts. With optically removable or controllable coatings, 
one could very conveniently control the functioning of the 
microdevice without acting directly on it. In addition, 
properly designed surfaces might provide temperature- 
controlled confining potentials which might find applica- 
tions in self-assembly processes [ooj . 
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